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Abstract 


The objective of this study is to develop a new methodology for estimating the reliability 
of engineering systems that encompass multiple disciplines. The methodology is 
formulated in the context of the NESSUS probabilistic structural analysis code, 
developed under the leadership of NASA Lewis Research Center. The NESSUS code 
has been successfully applied to the reliability estimation of a variety of structural 
engineering systems. This study examines whether the features of NESSUS could be 
used to investigate the reliability of systems in other disciplines such as heat transfer, 
fluid mechanics, electrical circuits etc., without considerable programming effort specific 
to each discipline. In this study, the mechanical equivalence between system behavior 
models in different disciplines are investigated to achieve this objective. A new 
methodology is presented for the analysis of heat transfer, fluid flow, and electrical 
circuit problems using the structural analysis routines within NESSUS, by utilizing the 
equivalence between the computational quantities in different disciplines. This technique 
is integrated with the fast probability integration and system reliability techniques within 
the NESSUS code, to successfully compute the system reliability of multi-disciplinary 
systems. Traditional as well as progressive failure analysis methods for system reliability 
estimation are demonstrated, through a numerical example of a heat exchanger system 
involving failure modes in structural, heat transfer and fluid flow disciplines. 
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CHAPTER I 


INTRODUCTION 

The purpose of this study is to develop a methodology to estimate the reliabil- 
ity of engineering systems that encompass several disciplines. The methodology is 
implemented using the NESSUS probabilistic analysis code, which has mostly been 
applied exclusively in the discipline of structural engineering. In order to apply 
the NESSUS probablistic structural analysis code to analyze a multi-disciplinary 
engineering system, the equivalences between system behavior models in different 
disciplines are investigated, and the effect of physical interaction among the failure 
modes is quantified in this study. 

System reliability analysis is a method of estimating the effects of uncertainties 
in an engineering system on the probability of successful performance. Usually, 
an engineering system consists of multiple subsystems and components, which may 
require the knowledge of different disciplines of engineering. Such disciplines may 
include structural engineering, mechanical engineering, heat transfer theory, fluid 
mechanics, electrical engineering, etc. Such a system is called a multi-disciplinary 
engineering system. The reliability analysis of any engineering system usually begins 
with the identification and reliability computation of individual failure modes within 
the system. Then the reliability analysis of the overall system can be carried out. 

Traditionally, relia’ iiity method?’ have primarily concentrated on failures in one 
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particular discipline, e.g. structural analysis, not on an overall system which con- 
sists of multiple disciplines. Furthermore, conventional methods of system reliability 
estimation usually only consider the statistical correlation between individual failure 
events, ignoring the fact that more often than not, those individual failure modes also 
have a physical correlation. This leads to inaccuracy in system reliability estimation. 

The method presented in this report uniquely computes the failure probability 
interactions between different modes and overal system failure probability through 
the imposition of one failure mode on another field and reanalysis of the latter. This 
method is used to compute the probabilities of critical system failure events after 
accounting for the contributing non-critical failure modes in all different fields. How- 
ever, it is not an easy task to estimate the reliability interactions between different 
failure modes. The success of such a method primarily depends on the availability of 
effective reliability tools. The software system NESSUS developed under the leader- 
ship of NASA Lewis Research Center is uniquely suited for this purpose. Currently, 
this code has been applied primarily to the structural engineering problems. In or- 
der to perform system reliability analysis including the interactive failure modes, this 
study uniquely develops behavior analogies between the structural model and heat 
transfer model, and between the structural model and fluid mechanics model. By 
doing so, the probability estimation of heat transfer and fluid mechanics failures can 
be pursued similarly to structural reliability analysis. 

The objective of this research project is to develop a method, using system relia- 
bility theory, for the reliability estimation of multi-disciplinary engineering systems. 
The method is implemented on the software ystem NESSUS (Numerical Evalu- 



ation of Stochastic Structures under Stress) developed by NASA Lewis Research 
Center. An example application to a three-discipline system involving mechanical 
stress-strain behavior, heat transfer and fluid mechanics is provided. In order to 
compute the individual failure mode probability of non-structural problems such as 
heat transfer and fluid mechanics within NESSUS, it is necessary to develop a new 
methodology for the analysis of heat transfer problems using the concept of equiva- 
lence between the computational quantities in structural analysis, such as stiffness, 
displacement vector, load vector, etc. and similar quantities such as conduction, 
temperature distribution and heat flux in heat transfer theory, flow velocity, pres- 
sure and flow factor in fluid mechanics. This is the first important contribution of 
this study. 

The second important contribution is the method for the computation of the 
physical dependence of critical failure mode probabilities on non-critical failure modes 
in various disciplines. This involves the imposition of the non-critical modes and 
reanalysis of the system with appropriate discipline equivalences, for various levels of 
progressive damage. The combination of these two ideas - inter-disciplinary analogies 
and physical failure mode correlation - makes a reliability analysis program such as 
NESSUS very powerful for application to a variety of multi-disciplinary systems. 

The concepts and methods discussed above are examined in detail in the next 
four chapters of this report. In Chapter II, the basic reliability analysis concepts 
for individual component-level and system-level events are reviewed, and their im- 
plementation in the NESSUS program is described. In Chapter III, the behavior 
analogies between the structural analysis model and heat transfer problem, and be- 



tween the structured analysis model and fluid mechanics model are developed. The 
finite element numerical examples with NESSUS/FEM are demonstrated for this con- 
cept. Chapter IV consists of two major parts: in the first part, the failure probability 
analyses for individual events including structural, heat transfer and fluid flow failure 
modes are performed using NESSUS; in the second part, the system failure proba- 
bility is studied. The effect of non-critical failure events of heat transfer and fluid 
mechanics upon a critical structural failure event is investigated, followed by system 
reliability analysis with the consideration of physically correlated component-level 
events. A numerical example of system reliability analysis of a multi-disciplinary 
system consisting of structural, heat transfer and fluid mechanical modes is demon- 
strated. The conclusions and recommendations of the study are summarized in 
Chapter V. 
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CHAPTER II 


SYSTEM RELIABILITY ANALYSIS 
Individual failure modes and effects 

An engineering system consists of a number of functional components. Before the 
system-level analysis begins, the modes of failure for individual components should 
be specified. The analyses of the failure modes and effects can be carried out by 
starting at the component level and expanding upward to the whole system. A 
failure mode is the manner by which a failure is observed. All units in a system are 
designed to fulfill one or more functions. A failure is thus defined as non-fulfillment of 
one of these functions. Analytically, each failure mode has a corresponding limit sate 
which separates the design space into “failure” and “safe” regions. The probability 
of failure, P/, is denoted as 

Pi = Pfe < 0] (1) 

where g is the value of the performance function g{X). The limit-state is denoted 
by the equation g(X) = 0. 

An exact solution of Pj can be obtained by the integration of the multiple integral 
denoted as 

Pj = / fx(x)dx (2) 

Jo(X)< 0 
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where f(X) is the joint probability density function of the vector of uncertain vari- 
ables X. 

In general, the solution of this multiple integral is too complicated to obtain. 
This is not only because the individual distributions axe not always available but 

T 

also because the integral is multi-dimensional for a realistic problem and is difficult 
to evaluate. Therefore, for practical purposes, efficient approximate analysis tools 
are needed. 

Fig. 1 illustrates the concept of the first-order approximation to the limit state 
for an estimate of the failure probability. 


The uncertain variables (X) are all transformed to equivalent uncorrelated stan- 
dard normal variables (u). The most probable point MPP of the limit state is defined 
at the minimum distance j3 from the origin to the limit state surface. Therefore, the 
first-order estimate of the failure probability is 

P. t = *(~0) (3) 


where $ is the distribution function of a standard normal variables. 

In the NESSUS computer code, this is referred to as the Fast Probability Inte- 
gration ( FPI ) method. The limit state is constructed as: 

g=Z(X)-Z o = 0 (4) 
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Figure 1: Failure Probability Estimation 


7 


where Zq is a real value of the random variable Z(X), which is a performance function 
or a response function, such as stress, displacement, temperature, etc. 

Z(X) = Z(X 1 ,X 2 ,...,X n ) (5) 

where X,(i = 1,2, ...,n) are the input random variables. 

The NESSUS program searches for the MPP by computing the sensitivities of the 
limit state to the random variables using iterative perturbation (in NESSUS/PFEM), 
and using these sensitivities to obtain a mean value first order ( MVFO ) or second- 
order {MV SO) estimate of the failure probability (in NESSUS/FPI). 

By using a first or second-order Taylor’s series expansion around the MPP , u*, 
the exact g(u)-function is replaced by the first-order polynomial, gi{u), 

n 

g x {u) = a 0 + ^2 cti(ui - u*) (6) 

«=i 

or a second-order polynomial, gziv), 

g 2 {u) = a 0 + f2 a»(«i - <) + 2 M u i “ u i ) 2 + S Y, Cij{ui - u*){uj - u*j) (7) 

i=l i=l 1=1 j= 1 

where the coefficients can be obtained by perturbation. 

Once these functions are obtained, the MPP is found. The probability of failure 
can be computed easily using Equation 3. This is the mean value first order {MV FO) 
estimate of the failure probability. This is improved using the Advanced Mean Value 
{AMV) analysis. Point probability estimate is made using specific limits for Z 0 , and 
the cumulative distribution function (CDF) is obtained by varying Zq. 
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Figure 2: Series System Structure 

NESSUS/FEM employs innovative finite element technology and solution strate- 
gies. It provides a choice of algorithms for the solution of static and dynamic prob- 
lems, both linear and nonlinear, together with an interactive perturbation analysis 
algorithm to evaluate the sensitivity of the response to small variations in one or 
more user-defined random parameters. 

NESSUS/FPI (Fast Probability Integrator) is used to evaluate structural response 
cumulative distribution functions (CDF). There are two methods in the code, the 
first-order reliability method and the advanced first-order reliability method. In gen- 
eral, the structured performance or response functions (e.g., stresses, displacements, 
vibration frequencies) are implicitly defined and each function evaluation may require 
intensive computation. The AMVFO ( Advanced Mean Value First Order) method 
reduces the computational burden and is the main probabilistic tool in NESSUS. 
NESSUS/PFEM automates the AMV procedure by integrating the FPI code and 
the FEM code. 
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System and component-level failure modes 


After the individual reliability analysis is completed, one then proceeds to the 
system or subsystem level analysis. System failure may occur due to a combination 
of any of the individual component failure modes. Many physical systems that are 
composed of multiple components can be classified as series-connected or parallel- 
connected systems, or combinations of series and parallel conditions. Description of 
these simple system structures is as follows. 

• Series System 

A system that is functioning if and only if all of its n components are functioning 
is called a series system structure. Fig. 2 illustrates such a system. 

If Ei denotes the failure mode i, then the failure of a series system is the event 

Ej = Ei U E 2 U ... U E n (8) 

Then the failure probability of the system is 
P, = P(E,) (9) 


If each failure mode Ei is represented by a limit state g(X) = 0 in basic variable 
space, the failure probability can be obtained by the integration denoted as 


= / ■•• / fx(x)dx 

J X J 
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(10) 




a 



b 


Figure 4: Parallel System Structure 

where X represents the vector of all the basic random variable (loads, material 
properties, etc.) and Cl is the domain in X defining failure of the system. 

This is defined in terms of the various failure modes as gi(X ) < 0. In two- 
dimensional X space, expression (10) is defined in Fig. 3. 

• Parallel System 

A system that is functioning if at least one of its n components is functioning is 
called a parallel system structure. A parallel structure of order n is illustrated in 
Fig. 4. 
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Figure 5: Combined System Structure 
In this case, the system failure event can be written 

Ej = Ei n e 2 n ... n E n (11) 

• Combined Series-Parallel System Structure 

This refers to systems which axe a combination of series and parallel structures. 
Fig. 5 shows an example of such systems. 

The failure event of this system is written, for example, as 

Ej = [Ei fl (E 2 U £ 3 )] U £4 ( 12 ) 

It should be noted that not all engineering systems cam be represented simply as 
described above. Practical systems may be more complex and need more effort to 
model. 
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System reliability computation 

In NESSUS, two methods have been implemented for system reliability compu- 
tation [4]: (1) probabilistic fault tree analysis combined with importance sampling 
(Torng et al, 1992) and (2) a structural reanalysis procedure to accurately estimate 
the failure regions for various critical failure modes affected by progressive damage 
(Mahadevan et al, 1992). 

Consider an engineering system subject to a sequence of loads (duty cycles) and 
which may fail in any one (or more) of a number of possible failure modes under any 
one load in the loading sequence. The total probability of the system failure may 
then be expressed in terms of the individual mode failure probabilities as 

Pj = P(Ei ) u p(e 2 n Si) u p(e 3 n s 2 n s x ) u p(e 4 n s 3 n s 2 n Si) u ... (13) 

where E{ denotes the “failure of the system due to failure in ith mode and S, denotes 
the complementary “survival event of the ith mode. 

Since P(E 2 fl Si) = P{E 2 ) — P(E 2 fl E\), Eq. 13 may be written also as 


Pj = P(Ei) + P(E 2 ) - p(Ei n e 2 ) + P(E 3 ) - P(Ei n e 3 ) 

-P{E 2 n E 3 ) + P{Ei n E 2 n e 3 ) + ... (14) 

where {E\ fl E 2 ) is the event that failure occurs in both modes 1 and 2, etc. 

Since it is not always an easy task to determine the joint probabilities of more 
than two failure modes, the following approximation methods can be used to predict 
the system reliabilities. 
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• First-order bounds 


The probability of failure for the system can be expressed as Pj = 1 — P(S), 
where P(S ) is the probability of survival. For independent failure modes, P(S) can 
be represented by the product of the mode survival probabilities, or, noting that 

P(Si ) = 1 - P(E { ), by 

Pj = 1 - fill - JW] (15) 

1=1 

where, as before, Pj is the probability of failure in mode z. This result can be shown 
to be identical with Eq. 14. It follows directly from Eq. 14 that, if P(Ei) 1, then 
Eq.15 can be approximated by [Freudenthal et al., 1966] 

P/ = E P(E<) (16) 

i=l 

In the case where all failure modes are fully dependent, it follows directly that the 
weakest failure mode will always govern system failure, irrespective of the random 
nature of the strength. Hence 

P, = max^S)] (17) 

Equations 15 or 16 and 17 can be used to define relatively crude bounds on 
the failure probability of any system of the series types when the failure modes are 
neither completely independent nor fully dependent. These are Cornell’s first-order 
bounds: 5 

max />(£,) < P(UUE.) <t,P(Fi) (IS) 
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• Second-order bounds 


For some practiced systems, the above first-order bounds may be too wide to 
be meaningful. For more accurate estimation, second-order bounds have been de- 
veloped. There are various second-order bounds in the literature [Kounias, 1968; 
Vanmarcke, 1973; Hunter, 1976; Ditlevsen, 1979]. Cruse et al (1992) derived second- 
order bounds which are independent of any ranking of the failure events[l]. The 
upper bound is 


^(U" =1 £,) < {£P(£,)-rr, M (£rnaxP(£,E J )max £ P(£, ■£,•)]} (19) 

t = l »=2 1=1, 

The lower bound is 


f(ur =1 £,)> m<K {/>(£,)+ J 2 ™x{(P(Ei)- £ 'p(EiE h )); 0 }} ( 20 ) 

k=l,k?i 

Utilization of the second-order bounds requires evaluation of terms of the form 
P(EiEj) where E{ denotes the event ’’failure in limit state The intersection terms 
refer to domains such as fix shown bounded by the non-linear limit state functions 
gi(X) = 0(i = 1,2,3) in Fig. 3. The individual failure mode probabilities in the 
first-order analysis are determined as 


P f = <K-0) (21) 

In standardized independent normal X space, the linear limit state function is 
given by 

71 

9i{X) = /?; + J2 °H x i (S2) 

j= i 
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where n is the number of random variables. 


The angle between the two limit states provides information about the correlation 
of the two failure modes. The correlation coefficient is obtained as 

n 

Pij = dirOjr = COSUij (23) 

r= 1 

Once Pj and are obtained, the computation of joint probability of failure 
can be carried out. Eg. 19 and 20 can be used to compute the second-order bounds 
for system failure probability estimation. 

The above method only provides a tool to approximately estimate the failure 
probability correlation of two different failure modes in a multi-discplinary system. A 
more accurate approach would be the imposition of one failure mode on anther mode 
and reanalysis of the latter. For example, consider two failure modes in a heat transfer 
system: structural failure and heat transfer failure. Structural failure happens when 
the stress, caused by the fluid pressure and temperature difference between outer and 
inner surfaces exceeds limiting value of strength. The heat transfer failure happens 
when the temperature of the contained liquid can not be kept at a certain level. 
When the thermal failure occurs, the increase in the temperature field also causes 
changes in stress field. The structural failure probabilty can be re-estimated under 
this changed stress field and the result can be considered. as the interactive failure 
probability under influence of heat transfer failure. A numerical examples will be 
shown for this approach in a later chapter. 
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Symbol 

Meaning of symbol 


Representation of an event 


Representation of an event of a failure 

Q 

OR-gate 

— 

AND-gate 


Figure 6: Standard symbols used in fault tree analysis 

Probabilistic fault tree analysis 

NESSUS system risk assessment (SRA) uses probabilistic fault tree analysis 
(PFTA). A fault tree is a mathematical construction of assumed component fail- 
ure modes (bottom events) linked in series or parallel leading to a top event, which 
denotes system failure. Standard graphical symbols are used to construct the fault 
tree picture, by describing events and logical connections. These axe shown in Fig. 
6, and a simple PFTA is shown in Fig. 7. 

• Fault Trees with a Single AND-GATE 

Consider the fault tree in Fig. 8. Here the top event occurs if and only if all 
the bottom events E\, E 2 , ..., E n occur simultaneously. A system with AND-GATE 
is very similar to a series system structure. 


18 



















Symbols 

AND-gate 


Description 



The AND-gate indicates that the output 
event A occurs only when all the input 
events Ej occur. 


OR-gate 



The OR-gate indicates that the output 
event A occurs if any of the input events 
Ej occur. 


Figure 8: Fault tree with a single AND-GATE and a single OR-gate 
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Figure 9: Schematic of NESSUS 

• Fault Tree with a Single OR-GATE 

Consider the fault tree in Fig. 8. The top event occurs if at least one of the 
bottom events E\, Ei , ..., E n occurs. The structure of this fault tree is similar to the 
paralell system structure. 

A schematic of Version 6.0 of the NESSUS (Numerical Evaluation of Stochas- 
tic Structures Under Stress) probabilistic structural analysis computer program is 
shown in Fig. 9. As shown in the diagram, the the NESSUS includes other modules, 
namely the System Risk Assessment (SRA) and Simulation Finite Element (SIM- 
FEM) modules. The random field pre-processor (PRE) provides data manipulation 
needed to express the uncertainties in a random field as a set of uncorrelated random 
variables. The user-subroutine which defines the response model (UZFUNC) enables 
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users to define required limit state with the computed response. This study will 
mainly use FEM, PFEM and FPI for reliability analysis. The NESSUS program is 
quite comprehesive with respect to structural reliability estimation. As mentioned in 
Chapter I, the purpose of this study is to develop a technique by which the NESSUS 
program can be used for the system reliability analysis of multi-discplinary systems. 
The following chapters describe this technique in detail. 
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CHAPTER III 


ANALOGY BETWEEN ENGINEERING SYSTEMS 

Introduction 

Since NESSUS/FEM program has been mostly applied only to structural analysis, 
a thermal or a fluid mechanical system needs to be converted through an analogous 
model to a structural system on which the NESSUS program cam be applied for 
analysis. Then the probability analysis for a heat transfer system or a fluid mechanics 
system can be carried out by NESSUS. By doing so, a system with heat transfer, 
fluid flow and mechanical stress problems can be analyzed by NESSUS automatically 
with FEM, PFEM, FPI and SRA modules for system reliability analysis. 

In this chapter, a new methodology is presented for one-dimensional steady-state 
heat transfer analysis and one-dimensional steady-state uniform flow problem using 
a structural finite element program. First, the use of the analogous models is intro- 
duced for the analysis of systems involving one-dimensional steady-state heat transfer 
and simple one-dimensional steady-state uniform flow in closed conduit systems. 

Heat transfer analysis through structural analogy 
One-dimensional steady-state heat transfer 

We begin our analysis of one-dimensional, steady-state conduction by discussing 
heat transfer with no internal generation. The objective is to determine the expres- 
sions for temperature distribution and heat transfer rate in common geometries. 
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Figure 10: Differential volume for the derivation of the general equation of heat 
conduction 

The concept of thermal conductivity (analogous to stiffness in stress analysis) is 
introduced as an aid to solving conduction heat transfer problems. Consider a three 
dimensional differential volume shown in Fig. 10, The general heat equation is 


d ,„dT. , d ,„dT. d ,„dT x . dT 
dx^ + dy^ K dy^ + dz^ K dz^ + q ~ pCp dt 


dx 


(24) 


where K is the thermal conductivity of the material. are related to 

heat flux in a direction perpendicular to the surface, q is the rate at which energy is 
generated per unit volume of the medium. The density p and specific heat c p are two 
thermodynamic properties. The product pc p is the volumetric heat capacity. pc p ^ 
is the time rate of change of the internal (thermal) energy of the medium per unit 
volume. 
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Figure 11: Heat transfer through a plane wall 

If the heat transfer is one-dimensional and steady state, any differentiation with 
respect to time is equal to zero and there is no internal heat generation, so Eq. 24 
reduces to 


d_ 

dx 




(25) 


The heat flux is a constant, independent of x. 

As shown in Fig. 11, a plane wall separates two fluids of different temperatures. 
Heat transfer occurs by convection from the hot fluid at T^i to one surface of the 
wall at T 4| i , by conduction through the wall, and by convection from the other surface 
of the wall at T, i2 to the cold fluid at T^. 

Assuming the thermal conductivity of the material to be constant, Eq. 25 may 
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be integrated twice to obtain the general solution 


T{x) = C x x + C 2 


To obtain the constants of integration, C x and C 2 , boundary conditions must be 
introduced. These are: 


T(0) = T 4i1 
T(L) = r , l2 


Applying the condition at x = 0 to the general solution, it follows that 


T s ,i = C 2 


Similarly, at x — L 


T,,2 = C X L + C 2 = C X L + T, iX 


in which case 


T a , 2 - Ts, 1 


= Ci 


Substituting into the general solution, the temperature distribution is then 

r <*>-K ^ ] { ?‘ } 

where N x = 1 — f, N 2 = £ 

The heat flow can be determined by Fourier’s law, that is 

dT 

"=- ka t x 


= -KA 


-(Ti - T t ) 

L 


— = - Ti) 
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Figure 12: Bar subject to tensile force F 
Stress analysis of a bar element 

Now consider a linear-elastic, constant cross-sectional area (prismatic) bar element 
shown in Fig. 12. Using Hooke’s law, the differential equation governing the linear- 
elastic bar behavior is 



= 0 


(34) 


where U is the axial displacement function in the x direction and S and E are 
cross-sectional area and Young’s modulus of elasticity respectively. 

M* *]{&} 

where Ni = 1 — and N 2 = 

The strain-displacement relationship is 

dU _U 2 - U x 
£x ~ dx~ D 
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Table 1: Analogous quantities for structural and thermal systems 


Heat transfer 

Structure 

Heat flux q 

Nodal force fi 

Temperature T(x) 

Displacement U(x) 

Inverse of heat transfer resistance 
Conduction: ^4, convection: hA 

Structural stiffness ^ 


We obtain 

F = ES (36) 

Also, by the nodal force sign convention of Fig. 12, 
h = ~F (37) 


So Eq. 36 becomes 
F <7 

h = ~ U 2 ) (38) 

Analogous modeling between heat transfer and structure 

Comparing Eq. 38 with Eq. 33, the similarities become apparent. These two 
equations indicate a direct analogy between heat transfer and structural analysis. 
The analogous quantities are listed in the Table 1. 

With this analogy, we are able to model a heat transfer problem into a stress 
analysis problem. 
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In the plane wall, we refer heat transfer resistance of conduction R to that 


is 


Rcond — 


t 3 , 1 - r Ji2 


KA 


(39) 


Considering the structural system, Hooke’s law provides stiffness of the form 


k= h = ES 

Ui-U2 D 

Comparing Eqns. 39 and 40, and considering ^ and ^ as analogous qualities, 
can be traced to be analogous to K. 

A heat transfer factor may also be associated with convection at a surface. From 
Newton’s law of cooling, 



g = hA(T, - T^) 


(41) 


where h is Planck’s constant of convection heat transfer coefficient, T, is the surface 
temperature and Too is the ambient temperature. 

The thermal resistance for convection is then 

D T, - Too 1 

ttconv — — * A 

q hA 

The equivalent thermal circuit for the plane wall with convection surface condi- 
tions is shown in Fig. 11. The heat transfer rate may be determined from seperate 
consideration of each element in the network, that is, 



q — hiA(Too,i — T, t i) 


KA 


(T aA - T, a ) = h 2 A(T 3<2 - Too, 2 ) 


( 43 ) 
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In terms of the overall temperature difference, To^i — t2 , and the effective 

thermal resistance R e //, the heat transfer rate may also be expressed as 
Too,l ~ Too, 2 


<? = 


R, 


*11 


(44) 


Because the conduction and convection resistance are in series and may be summed 
up, it follows that 


D 111 
R*/f = 7— r + "7TT + 


(45) 


h x A KA h 2 A 

Consider a bar consisting of three different materials which are denoted as ele- 
ments 1,2 and 3. The effective stiffness for this composite bar is 


k *ll = XT TIT 

ki ^ fc 2 T fe 3 


(46) 


Comparing the above equations Eq. 45 and Eq. 46, the analogy is k e jf r R ^ , 
that is 


111 
ReJf < ~^r l + r 2 + k^ 

Substituting with Eq. 45, we obtain 


(47) 


1 L 1 

+ T7~Z + 


hxA KA h 2 A 
where 

Ei Si 


J_ J_ J_ 
h + k 2 + k 3 


(48) 


k x = 


D x 


(49) 


k 2 = 


KiSx 

Do 


(50) 


k* = 


EsSx 

Do 


(51) 
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Substituting the corresponding terms in Eq.48, we obtain the equivalent quanti- 


ties 


E i 



(52) 


E 2 



(53) 


E 3 (V*)(|r) (54) 

With these analogous quantities, we use the NESSUS/FEM beam element with 
the E values replaced by the values involving heat transfer problem outlined above. 
The boundary conditions for the bar are the end displacements corresponding to the 
ambient temperature of the wall. After the structural analysis, we get the tempera- 
ture distribution from the corresponding displacement distribution in the output. 

Heat transfer in composite walls 

Equivalence concepts for thermal-structural analysis may also be used for more 
complex systems, such as composite walls and radial heat transfer systems. Fig. 13 
shows a series composite wall. The one-dimensional heat transfer rate for this system 
can be expressed as 


<7 = 


Too ,! ~ Top, 4 

T.R 


(55) 


where T 00i i — is the overall temperature difference and the summation includes 

all thermal resistances. Hence, 
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I I I T 00 , 4 , h 4 

Cold fluid 


h,A k A A k B A k c A h 4 A 

Too, 1 Ts, 1 T2 T3 Ts ,4 Ao °’ 4 


Figure 13: Equivalent thermal circuit of a series composite wall 
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A 


Figure 14: Structural analog for the series composite wall heat transfer 


<1 = 


Toq,i — TqQ'4 

Ri + R? + R 3 + R 4 + Rs 


(56) 


7 qq,i — Too t 4 

(1 /hx A) + ( L a IK a A ) + ( L b IK b A ) + ( L B /K B A ) + (l/h 4 A) 


(57) 


Alternatively, the heat transfer rate can be related to the temperature difference 
and resistance associated with each element. For example, 


T i rri m rp rrt rrt 

00,1 1 9 t 1—^2 _ 2 ^3 

(1/AiA) “ (L a /K a A)~ (L b /K b A) 


(58) 


The analogous structured model for this series composite wall heat transfer prob- 
lem is shown in Fig. 14. The bar consists of five elements with stiffnesses of 
fci, fc 2 , A: 3 , k 4 ,k 5 . Using the mechanical structure equivalence for convection and con- 
duction, we obtain 



(hiA) (§) 
(¥) (§) 


for convection 

► 

for conduction 
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Engine oil 



A- 


y?t Q) 2\3 (j) 4\s Q <1 6 


A 


(b) Structural model 

Figure 15: A heat exchanger for engine oil and refrigerant fluid 



Heat transfer in radial systems 


Cylindrical and spherical systems often experience temperature gradients in the 
radial direction only and may therefore be treated as one dimensional. 

Fig. 15 shows an example of a heat exchanger, whose inner cylinder is used to 
store engine oil and the outer cylinder is used to transfer the refrigerant fluid to 
cool down the oil temperature. The outer insulated covering is assumed to isolate 
the system from the ambient environment. For steady state conditions with no heat 
generation, the appropriate form of the heat equation is 



The rate at which energy is conducted across any cylindrical surface in the solid 
may be expressed as 


q = -KA— = -K(2nrL)^~ 
dr dr 


where A = 2nrL is the area normal to the direction of heat transfer. 
The thermal resistance is 


R 1 lnfa/rQ 1 

e// MTrrxL 2 nKL h 2 2nr 2 L 


(60) 


(61) 


which includes both conduction and convection. 

The heat transfer rate for a unit length of the cylinder therefore is 

Top,! — Tpo, 2 

^ _ 1 I ln ( r 2/ r l) . 1 

hi 2 nriL 2 nKL ' /12271T2L 


( 62 ) 
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Table 2: Analogous quantities for heat transfer in a radial system. 


Structure 

Heat transfer 

Ei 

(2*rihi) (§f) 

Ei 



E 3 

(2nr 2 hi) ( 7 ^) 


The structured analog for the cylinder is shown in Fig. 15(b). In this case, the 
mechanical structure equivalence is k e jj < — > The corresponding equivalent 

quantities are listed in Table 2. The E values are input to NESSUS structured anal- 
ysis, and the output displacements from NESSUS give the temperature distribution. 

Numerical example for heat transfer solved with NESSUS/FEM 

Fig. 17 shows the sectional view of the cylindriced copper heat exchanger which 
the engine oil flows through. The copper wall thickness is 0.281 in. The radius 
to the surface of the insulation pipe covering ( ki = 0.428 Btu/(h — ft 2 F )) is 

1.33 in. The fluid in the outer container is controlled at a constant temperature 
of 70 °F. The forced convection heat transfer occurs between the outer surface of 
the insulation covering and the flowing fluid with h = 10 Btu/(h — ft 2 — ° F). The 
surface temperature at the insulation covering is 35.298° F. The structural analogy 
model is used to determine the inside temperature of the tube, assuming steady 
state, one dimensional, uniform properties in each material, forced convection cooling 
ar'd negligible thermal a ..ation. The conductivity coefficient of copper at room 
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(a) Thermal circuit 


Jc3 = EjS j/D 3 ki = E. S i/Dj 



(b) Structural analogous model 


Figure 16: Analogous model for the heat exchange 

temperature is K c = 223 Btu/(h — ft 2 F). 

The thermal circuit is shown in Fig. 16(a). Fig. 16(b) shows the structural 
analog model for this heat transfer problem. Beam element type 98 in NESSUS/FEM 
element type library is adopted. Three elements represent three heat transfer forms 
involved in this problem, which are the forced convection between the surface of 
the insulation covering and the ambient air, conduction through the copper layer, 
and conduction through insulation covering, respectively. Therefore, in terms of the 
structural model, we must assign three different material elastic constants for this 
beam structure. Since the NESSUS/FEM utilizes the Nodal-based data input , two 
duplicate nodes are used at each boundary between elements 1 and 2, and between 
elements 2 and 3. The room temperature 70 °F becomes the boundary displacement 
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r 0 = 1.330 in 
Tj - 0.669 in 
r 2 = 0.950 in 


7. = 7(fF 

H- 10 Btu /(h - ft 2 -°F) 
k = 0.428 Btu/(h- ft 2 -°F) 


Figure 17: Section of a radial heat exchanger 


70.0 at node 1. A concentrated load F at point 3 is -241.66 lb. The length of the 
structural element is 10.0 units, the sectional area is 1.0. It should be noted that the 
units used in the structural model here do not have any read meaning in terms of a 
real structure. They are simply used to facilitate the structural analysis. 

The equivalent values are calculated as follows. 

( 2ttKc \ /£i\ 2 tt x 223 


El (.Si) ln(0. 95/0. 669) 


x 10 = 39955.476 


E, 


_ / 2 nKi 


Di 


) = 


2tt x 0.428 


\ln(r 0 /r 2 )/ \ 5 2 / ln(1.33/0.95) x 10 


E 3 = (2tt r 0 h) = 2tt x 11.33/12 x 10.0 x 10 = 69. 


= 79.923 


63867 


(63) 

(64) 

(65) 


F = 27rr 0 fc(r oo ~ T 3 ) = 2n x 1.33/12 x 10(35.298 - 70) = -241.66 


( 66 ) 
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NESSUS/FEM uses this data, and gives the output of the displacement distribu- 
tion in the structure as follows: 

U x = 70.000 
U 2 = 35.298 
U z = 35.298 
U A = 5.0614 
U 5 = 5.0614 
Ue = 5.0009 

Converting the above displacement information to the equivalent temperature 
distribution, we obtain: 

Too = 70.000° F 
Tx = 35. 298° F 
T 2 = 5.0614° F 
T 3 = 5.0009° F 

The data and the output files are shown in Appendix A. 

Fluid flow analysis through structural analogy 

Equation of motion for fluid flow 

The Bernoulli equation gives a relationship between pressure, velocity, and posi- 
tion or elevation in a flow field. Normally, these properties vary considerably in the 
flow, and the relationship between them if written in differential form is quite com- 
plex. The equation can be solved exactly only under very special conditions. There- 


39 



Figure 18: Fluid in a constant diameter duct 
fore, in most practical problems, it is often more convenient to make assumptions to 
simplify the descriptive equations. The Bernoulli equation for steady, incompressible 
flow along a streamline with no friction (no viscous effects) is written as [10] 

P V 2 

p + ~2~ +gz = C ( 67 ) 

where 

p is fluid pressure 
p is the desity of the fluid 
V is the flow velocity 
g = 32.174 ft/s 2 , and z is height. 

For a horizontal pipe shown in Fig. 18, zi = z 2 . From continuity, A x \ \ = A 2 V 2 . 
Because D i = D 2 , then A\ = A 2 , and therefore, V] = V 2 . The Bernoulli’s equation 
reduces to 

Pi = P'2 (68) 
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X w Pdz 



(p+dp)A 


-e- 


Figure 19: Control volume of a system: flow in a duct. 

This result is not a proper description of the situation, however. For flow to be 
maintained in the direction indicated in Fig. 18, p\ must be greater them p? in an 
amount sufficient to overcome friction between the fluid and the pipe wall. In order 
to apply Bernoulli’s equation and obtain an accurate description, we must modify 
the equation with a friction term. 

Consider flow in a pipe as shown in Fig. 19. A control volume that extends to 
the wall (where the friction force acts) is selected for analysis. 

Note that a circular cross section is illustrated, but the results are general until 
we substitute specific equations for the geometry of the cross section. The forces 
acting on the control volume are pressure normal to the surface and shear stress 
acting at the wall. The momentum equation is [10] 

£ F* = / J V z pV n dA (69) 
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where 

V z is the fluid velocity along the longitudinl direction 
V n is the norminal fluid velocity 


Since the flow out of the control volume equals the flow in, the right-hand side 
of this equation is zero. The sum of the forces is 

pA - r w Pdz — (p + dp)A = 0 (70) 

where 

A = cross-sectioned area 

Pdz = the surface area (perimeter times length) over which the 
wall shear t w acts 

The equation reduces to 


r w Pdz + Adp = 0 


(71) 


Rearranging and solving for pressure drop, we get 

dp _ _4 t w 
dz D h 


(72) 


We have thus expressed the pressure drop per unit length of the conduit in terms 
of the wall shear and the hydraulic diameter. Eq. 72 is a general expression for any 
cross section. It is convenient to introduce a friction factor /, which is customarily 
defined as the ratio of friction forces to inertia forces: 


/ = 


4t 

I > w 

W 2 


(73) 
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Figure 20: Laminar flow in an annulus 

where V is the average flow velocity. 

By subsitution into Eq. 72, we obtain 

J pV J Jdz 

dp= -—-D^ 

Integrating this expression from point 1 to point 2 a distance L apart in the 
conduit yields 



V 2 = p (75) 

pfL 

Eq. 75 gives the relationship between the velocity and the pressure drop in the 
duct due to friction. This equation can be applied to two flow regimes - laminar and 
turbulent flow. However, caution must be excercised when determining the friction 
factor /. 

This equation can also be applied to flow through noncircular cross section such 
as rectangular duct and annulus. Fig. 20 shows the laminar flow in an annulus. 
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Table 3: Antilogous quantities between structural and flow systems 


Fluid mechanics 

Structure 

Square of velocity V 2 

Nodal force fi = —F 

Pressure distribution p(x) 

Displacement U[x) 

Flow factor jj£ 

Structural stiffness ^ 


The annulus flow area is bounded by the inside surface of the outer duct (radius 
Ri) and the outside surface of the inner duct (Ri). We define the ratio of these 
diameters as 


k = 


R2 

Ri 


(76) 


in which 0 < k < 1. 

The friction factor used in Eq. 75 is defined as [10] 


l_ Re [ 1 ±Jf_ 1 +fc 
/ 64 1 — k + In (k) 


(77) 


where R e = e¥S2 E l ( 1 — fc). 

Compare Eq. 75 with Eq. 38 concerning the beam structure subjected to the 
end nodal force, as discussed in previous section: 


pc 

h = ^( U ! - U 2 ) 


(78) 


We are now able to set up the analogous quatities listed in Table 3. 
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Numerical example of flow in a tube solved with NESSUS/FEM 


Consider the refrigerant flow in a copper tube as an example to demonstrate how 
NESSUS/FEM can be applied to problems in fluid mechanics. 

A horizontal copper duct as shown in Fig. 18 with inside radius of 0.669 in, and 
1,200 in in length. If the inflow pressure p L is 1838.7 psi, assuming the refrigerant 
is Freon F-12 under a temperature of 5 °F, p is 0.0499 lb/ in 3 . The friction factor / 
is assumed to be 0.03, and V is 15.5 in/ sec. The objective is to obtain the outflow 
pressure p? using NESSUS/FEM. 

First of all, we need to identify all the equivalent quantities for structural analysis. 
We assume a single element beam structure subjected to a concentrated force equal 
to -240.25 units. The beam element has a section of 0.1 inx 0.1, and a length of 1.0. 
The boundary condition is an initial displacement of 1838.3 units at node 1. Again, 
it should be mentioned that the units used here do not have real meaning in terms 
of a real structure. According to Table 3, the analogous quantities can be obtained 
as follows 



2 x 2 x 0.669 
0.0499 x 0.03 x 1200 


x 10 = 14.896 


f x = V 2 = 15.5 2 = 240.25 


F = -240.25 

This data is input to NESSUS/FEM, and the displacement at point 2 is ob- 
tained as 1677.4. Converting this displacement to the fluid model, we get the output 
pressure p i — 1677.4 psi. 

The data and the output files are shown in Appendix B. 
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CHAPTER IV 


MULTI-DISCIPLINARY SYSTEM RELIABILITY ANALYSIS 

Introduction 

After transforming the heat transfer and fluid mechanics problems into corre- 
sponding structural analog models and using NESSUS/FEM to perform the finite 
element analysis, we can define the individual failure modes in NESSUS/FPI. Then 
NESSUS/PFEM can be employed to integrate FEM and FPI programs to obtain 
the failure probability and CDF for each failure mode. The failure mode for heat 
transfer problem would be defined as, for example, the event that the temperature 
at a certain location is lower or higher than the required temperature. The failure 
mode for fluid flow would be defined as the flow pressure exceeding a certain pressure 
level, and the structural failure is defined as the stress exceeding either the ultimate 
strength or the yield strength of the material. 

Upon the completion of failure probability analyses of individual failure modes, 
the system failure analysis can be pursued. The different failure modes involved in 
a system have different impacts on the overall performance of a system. Some types 
of failure such as structural failure are critical to the system. If the material used to 
construct the main parts of the system fails, the whole system can no longer function. 
Such failure is called critical failure. Other failures modes such as thermal failure of a 
hect exchager do not destroy the system but degrade the performance of the system. 
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Such failure is referred as to functional failure. The function of fluid flow will fail 


when the outflow pressure rises higher than the designed value, but the system can 
still be working until the pressure increases to the level which will cause the system 
to shut down. The individual failure modes can also be correlated to each other. For 
example, the temperature field in the thermal failure mode affects the stress field in 
the structural mode. The flow pressure definitely has impact on the stress. However, 
in some cases, the component-level events in a system is considered as independent 
events. In the example which will be discussed later in this chapter, the thermal 
failure and fluid flow failure modes do not share correlated input parameters, so they 
are considered as independent of each other. 

Using the analogy method, the thermal and fluid flow problems axe analyzed sim- 
ilar to the structural model by means of NESSUS. For physically correlated events, 
the failure mode of one event is imposed on the other. In a system consisting of 
structural, thermal and fluid flow modes, the thermal and flow failures are imposed 
into the structural failure analysis to study the impact of correlated events. The 
failure probability of the whole system is then estimated based on the output from 
the above analyses. 

Individual failure analysis 

Structural failure mode 

For a copper duct of a heat exchanger shown in Fig. 17 in Chapter III, the 
structural failure mode is defined as that when the tensile stress exceeds the yield 
strength f y i e id = 8.0 ksi. In finite element modeling, we use Element type 153 in the 
NESSUS/FEM file. This element is a four-noded quadrilateral lying in the global 
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zt— plane which is defined by cylindrical coordinates. 

This structure is subjected to two types of load - fluid pressure from the inside 
flow and the stress caused by temperature difference between the outer surface and 
the inner surface. 

One convenient feature of the NESSUS/PFEM is that we can impose different 
temperatures at the inner surface of the pipe and obtain the different probability 
results under different temperature conditions. This enables us to investigate the 
effect of different temperature levels on the structured failure probability. By doing 
so, the relationship between the failure modes in two dicsiplines - heat transfer and 
structural mechanics - is established. This is a significant step toward the system re- 
liability analysis with physically correlated failure modes. This will be demonstrated 
in a later section. 

For this structural failure model, we first suppose that the temperature failure 
(which will be descibed in the following section) did not occur, that is, the tempera- 
ture at the inner surface of the duct is below 5.0009°F. Given that the inner surface 
temperature is 4.0000°F, using the FEM we obtain the outer surface temperature as 
4.0605°F. 

Also, we assume that the outlet flow pressure is under 1677.4 psi which enables 
the system to work properly. We assume that the outlet flow pressure is 1577.4 psi. 

We input this temperature and flow pressure profile in the structure FEM data 
file, with the random variables defined in Table 4. 
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Table 4: Random variables for structural model 


Variable 

Mean value 

Distribution 

C.O.V. 

Pressure 

P 

1577.4 psi 

Normal 

0.2 

Modulus of elasticity 
E 


Normal 

0.1 

Coefficient of thermal 
expansion a 

9.5 xlO~ 6 

Normal 

0.1 


Table 5: CDF corresponding to different tensile strength levels 


Z-level 

(strength) 

CDF 

Z-level 

(strength) 

CDF 

-206.89941 psi 

0.00000017 

4385.4964 psi 

0.81593991 

558.49989 psi 

0.00002067 

5150.8957 psi 

0.97128351 

1323.8992 psi 

0.00096767 

5916.2950 psi 

0.99813412 

2089.2985 psi 

0.01786435 

6681.6943 psi 

0.99995188 

2854.6978 psi 

0.13566610 

7447.0935 psi 

0.99999952 
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The result attached in Appendix C indicates that the structural reliability when 
the heat exchanger is working properly in thermal and fluid aspects is 0.99999999. 
The failure probability is expressed as 1 — Preliability Therefore, the structural failure 
probability is 1.0 x 10 -9 . 

The key word response type in the NESSUS input data file, FPI section, is set 
equal to 3 which means that the response quantity used in limit state function is 
stress. The corresponding keyword analysis type in FPI section is first set equal to 
1 which means that the probability analysis is for a single Z-level. The Z-level in this 
case is 8,000 psi. The probability result will be under the condition of a < 8, 000 
psi, i.e., the structural reliability of the system under certain thermal and fluid flow 
working conditions. 

The CDF is obtained by using PFEM by setting analysis type in FPI section 
equal to 0 which automatically generates a set of different values of Zq (i.e., Z levels 
for a series of stress valus) for probability analysis. The CDF values corresponding 
to different strength Z-levels are shown in Table 5. 

The CDF chart is shown in Fig. 21. It should be noted that the first line of 
the data which contains negative Z-level is eliminated because negative stress is 
considered impractical in this model. The input and output files are attached in 
Appendix C as well. 
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CDF 


1.0 



Figure 21: CDF of structural reliability of refrigerant duct 


Table 6: Random variables for thermal model 



Mean value 

Distribution 

Coefficient of variation 

K c 

223.0 

Normal 

0.1 

Ki 

0.428 

Normal 

0.1 

h 

11.3 __ 

Normal 

0.1 
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Figure 22: CDF of internal fluid temperature of refrigerant duct 
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Failure mode in heat transfer 


Using the same example of a heat exchanger as in Fig. 17, we define a failure 
event when the inside temperature is higher than 5 °F, because the refrigerant will 
not function properly beyond 5 °F which is considered as a failure in the device we 
studied. First the input data for NESSUS/PFEM is set up to obtain the reliability 
under this failure mode, then the data file is set up with different Z - levels to obtain 
the CDF, which provides reliability estimatecorresponding to different temperature 
levels. The random variables K c and Ki and h for heat transfer are defined in Table 
6 . 

In order to use NESSUS/PFEM, the analogous quantities E\, F 2 , E$ and F 
are calculated from Eqs. 63, 64, 65, 66. Because the distribution of the random 
variables K c and K{ and h is normal and Fi, E 2 , E$ are linear to K c and Ki and 
h, the distribution of random variables E\, E 2 , E 3 is also normal. The mean values 
and standard deviations of of E\, E 2 and E 3 are input to NESSUS/PFEM. 

The Z-level is 5.0°F, so P(Z < Z 0 ) is the probability the device can keep the 
inside fluid temperature under 5.0°F, which is the thermal reliability of the system. 
We set up the keyword in FPI section analysis type equal to 1 which means the 
probability analysis is performed for a single Z-level. The result is attached in Ap- 
pendix D. The thermal reliability of this device is 0.9099214. Therefore, the thermal 
failure probability is 9.00786 x 10 -2 . 

The CDF is obtained by setting up the FPI keyword analysis type equal to 0 
which automatically generates a set of different values for Zq (i.e., Z-levels). The 
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Table 7: Random variables for flow model 



Mean value 

Distribution 

Coefficient of variation 

D h 

0.669 

Normal 

0.1 

V 

18.71 

Normal 

0.1 


CDF output is shown in Appendix D and the CDF curve is shown in Fig. 22. 

Failure mode in fluid flow 

Next we consider the one-dimensional fluid flow in a duct of a heat exchanger. 
The failure mode is defined as the pressure at a certain point along the duct rising 
above the value at which the system cannot function properly. 

The example of a duct in a heat exchanger shown in Fig. 16 is used. The only 
difference is that V is assumed to be 18.7 in/ sec. We define the failure mode when 
pressure rises above 1677.4 psi. The Z-level is therefore 1677.4 psi. The keyword 
response type is set as 1 for the displacement output which is the analogy of the 
pressure. The random variables related to fluid flow are defined in Table 7. The 
analogous quantities for use in NESSUS/PFEM are calculated according to Table 3 
as follows: 


E = 




2 x 2 x 0.669 
0.0499 x 0.03 x 1200 


x 10 = 14.896 


f l = V 2 = 18.7 2 = 350.0 


F = —fi = -350.0 
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Figure 23: CDF of fluid pressure of refrigerant in the duct 
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Since £ is a linear fuction of Dh, its distribution is also normal. But since f l -=V 2 1 
the distribution of f\ is actually chi-square (x 2 )- However, we have used the normal 
distribution for f\ in this study as an approximation. The friction fator / is assume 
to be a constant. 

A NESSUS/PFEM input data file is compiled. The reliability is obtained by 
setting the keyword analysis type equal to 1, and the CDF is obtained by setting it 
to a value of 0 which automatically generates a set of different values for Zq. The 
PFEM input and output files are shown in Appendix E and the CDF curve is shown 
in Fig. 23. 

The reliability is 0.98670241 and therefore the failure probability of the output 
flow pressure being higher than 1677.4 psi is 1.329759 x 10 -2 . 

Multi-disciplinary system reliability 

After the individual failure modes are identified and analyzed, the system reli- 
ability analysis can be pursued. Fig. 24 shows a device which is used to transfer 
refrigerant fluid through a copper duct. The duct is installed in an enclosed cham- 
ber which is maintained at a constant temperature of 70 °F. The thickness of the 
copper wall is 0.281 in. The radius to the surface of the insulation pipe covering ( 
mean value of K{ equals to 0.428 Btu/(h — ft 2 — 0 F), c.o.v. equals to 0.1 ) is 1.33 
in. Forced convection heat transfer occurs with h = 10 Btujih — ft 2 — 0 F) ( mean 
value with c.o.v. equals to 0.1). The thermal conductivity of copper K c has a mean 
value of 233.0 Btu/(h — ft 2 — 0 F) with a c.o.v. 0.1. The surface temperature at 
the insulation covering is 35.298° F. The inflow pressure p\ is designed to be 1838.7 
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psi. The refrigerant is Dichlorodifluoromethane ( Freon F-12 ) under a temperature 
of 4 °F. p is 0.0499 Ib/in 3 . The friction factor / is assumed to be 0.03, and V is 
18.71 inf sec ( assuming V 2 is normal distribution and has a c.o.v. of 0.1). The inner 
radius of the duct has a mean value of 0.669 in and a c.o.v. of 0.1. 

The above data axe the same as the model shown in Fig. 17 of Chapter III, which 
are used in FEM analysis and reliability estimation of individual failure modes. The 
system shown in Fig. 24 is simply the combination of the previous individual models 
which have been analyzed in different disciplines. The following is to demonstrate 
how the analysis results of the individual failure modes cam be integrated into the 
analysis of a whole system. 

The failure of the system consists of the individual failure modes in three disci- 
plines: structural failure, thermal failure and fluid mechanics failure. 

• First of all, the duct should work without any damage to the structure, i.e. the 
duct should be structurally sound without yied or crack. If yield occurs, then 
structural failure is assumed to occur. We denote the structural failure as E\. 
The structural failure is a critical failure in this system. 

• The refrigerant liquid this device transports is sensitive to temperature changes. 
The requirement is that the temperature cannot be higher than 5°F for the 
next process to proceed. If the temperature of the liquid rises higher than 5 °F, 
then thermal failure occurs, which we refer to as E 2 . The thermal failure is a 
non-critical functional failure in this system. 

a It is required that the fluid flow be maintained at a certain presure at the 
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Input valve 

Outlet to the next process 



Figure 24: Refrigerant duct through a chamber 
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ends which enables the refrigerant fluid to maintain a steady speed to provide 
constan volume in the device. If the flow pressure rises higher than 1677.4 
psi, the flow failure occurs, which is referred to as E 3 . The flow failure is a 
non-critical functional failure in this system. 

From the previous section, the probability of the individual failure modes of the 
system shown in Fig. 24 have been obtained, wich are: 

P{Ei) = 1.0000000 x 10 -9 
P(E 2 ) = 9.0078600 x 10’ 2 
P{E 3 ) = 1.3297590 x 10 -2 

As was indicated in the previous section, we can also impose one failure mode 
upon the other. In this case we can impose the thermal failure ( which happens 
when the fluid temperature' rises above 5 °F ) and the fluid pressure failure ( which 
happens when the outlet fluid pressure rises above 1677.4 psi) upon the structure 
respectively. In the FEM file for the structure model, the corresponding data are 
modified to impose those failures. 

First, we assume that the temperature failure occurs while the fluid pressure is 
still lower than 1677.4 psi , say 1577.4 psi, i.e., the fluid flow is operation in safe mode. 
By redefining the temperature profile in the FEM data deck as 6.06 ° F (for example) 
in the inner layer of the wall and 6.12°F in the outer layer of the wall, which mean the 
thermal failure occurs, we inpose the thermal failure to the structural model. The 
PFEM result gives us the structural failure probability under the condition that the 
thermal failure occurs. In this case, the structural reliability is 0.99999996, therefore. 
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the failure probability P{E\j E 2 ) = 4.0 x 10~ 8 . 

Now we impose the fluid pressure failure upon the structured model. The fluid 
mechanical failure occurs when the pressure at the outlet rises above 1677.4 psi, say 
1777.4 psi. The structural FEM file is modified by redefining the pressure profile 
according to this failure pressure. It should be noted that the temperature profile 
should remain under the normal working condition, which is that the temperature 
in the inner layer of the wall of duct is under 5°F, say 4°F, i.e., the thermal as- 
pect of the system is operating in the safe zone. The result indicates that in this 
case, the structural reliability is 0.99999633, therefore, the conditional probability, 
P(Ei/Es) = 3.67 x lO" 6 . 

Next, both thermal and fluid mechanical failures are imposed that is, the fluid 
temperature rises above 5 °F, and the outlet flow pressure rises above 1677.4 psi. 
Modifying the input FEM data deck in structural PFEM file with inner surface 
temperature of 6.06 °F, and the fluid pressure of 1777.4 psi, we can get the result of 
the structural reliability of 0.99999227, which means, P(E\/ E 2 E 2 ) = 7.73 x 10~ 6 . 

The conditional probabilities of structural failure have been obtained as 

P(E l /E 2 ) = 4.00 x 10‘ 8 

P(Ei/Ez) = 3.67 x 10~ 6 
P{E 1 /E 2 E 3 ) = 7.73 x 10" 6 

System reliability computation 

System reliability analysis can be performed in two different ways, depending on 
the definition of the systen failure. In the first (traditional) method, we define that 
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Structural Thermal Fluid flow 
Failure Failure Failure 


Figure 25: Fault tree for the system with three critical failures 
the system failure occurs when any component-level failure occurs. In the system 
invoving structural, thermal, fluid mechanical failure modes, i.e. , E\, E 2 and £ 3 , the 
system failure can be illistrated in a fault tree shown in Fig. 25. 

As discussed in Chapter II, the probability of system failure P(E) can be obtained 
using the following equations: 

£(£) = P(E X U £ 2 U E 3 ) (79) 

The above expression can be expaned as: 

P(E) = P{Ei) + P(E 2 ) - P(E \ n E 2 ) + P(E 3 ) - P(E X n E 3 ) 

— p(E 2 n £ 3 ) -h £(£i n E 2 n £ 3 ) (so) 

Since the joint probability is not always available, an approximate method is to 
consider the individual failure modes as independent and ignore the correlations. In 
our case however, the conditional prob:- bilities have been calculated. Therefore, the 
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system failure probability can be computed as 


P(E) = P(E X ) + P{E 2 ) + P(E 3 ) - P(E 1 /E 2 )P(E 2 ) - P(E l /E 3 )P(E 3 ) 

-P(E 2 /E 3 )P(E 3 ) + P(E l /E 2 E 3 )P(E 2 E 3 ) (81) 

For systems involving many failure modes, approximation methods are used to 
predict the system failure probability (or reliability), such as first-order bounds or 
second-order bounds [1]. 

Since no correlation is assumed between the failures of flow pressure and fluid 
temperature, we assume P(E 2 fl E 3 ) is equal to P(E 2 )P(E 3 ) in our analysis. 

Substituting the numerical results from the previous discussion into Eq. 81, 
we obtain the probability of the system failure of the heat exchanger, P(E) = 
0.10217832. 

As mentioned before, the above failure probability is an estimation of system 
failure in case any failure occures which includes both critical and non-critical func- 
tional failures. Now we will pursue the probability estimation for the system critical 
failure which, in our case, is structural failure. During the service cycles, the thermal 
and fluid mechanical failures may be non-critical, i.e., their occurence does not cause 
total system failure. They will cause the system to fail in some functions as designed, 
such as keeping the fluid under certain temperature or keeping the outlet fluid pres- 
sure under certain value. However, if the the system keeps operating, the changes in 
temperature and fluid pressure will cause progressive damage to the structure due 
to load redistribution. The estimation of critical structural failure of the system has 
to consider the progressiva damage caused by all components. Structural reanaiysis 
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Table 8: Structural failure probability under various temperatures 


TO 

0.0 

5.0 

10.0 

15.0 

P(E i) 

0.0 



0.0 


is used to account for the effect of non-critical damage on critical failure mode. In 
the refrigerant model, a reanalysis procedure is performed to accurately estimate the 
failure region segments for structural failure mode affected by progressive damage 
caused by thermal and fluid pressure changes within the system. The overall struc- 
tural failure probability is obtained through the union of the failure region segments 
defined by each limit-state function. 

We can also impose a series of temperatures under which the system may be 
operating upon the structural model to examine the temperature impact on the 
structural failure probability. Just as we did before, the failure probability is obtained 
as (1 - Reliability). In this case, we still assume inner fluid pressure is 1577.4 psi, 
which means that the fluid flow mode of the system is operating in the safe zone. 
The results are shown in Table 8. 

Table 8 shows that when the fluid pressure is not considered as a random variable 
in perturbation for probability analysis, the temperature changes do not have a 
significant impact on structural reliability of the system. 

We can also get the structural failure probability under different pressure condi- 
tions by defining a series of the pressure profiles in the FEM data deck for strucutural 
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Table 9: Structural failure probability under various pressures 


Pressure 

(psi) 

2700.0 

2900.0 

3000.0 

3200.0 



3420.0 

P(E X ) 

0.000 

1.000 

xlO- 9 

1.920 

xlO -6 

7.687 

xlO -3 

9.195684 

xlO" 2 

4.06468 

xlO -1 

4.92622 

xlO" 1 


model. The results are listed on Table 9. In this case, we assume that the inner tem- 
perature is 4°F. 

It should be noted that the NESSUS/PFEM input file gfun.dat is different from 
the previous structural PFEM file in which the pressure is defined as a random 
variable. In the program [Ang and Tang, 1984] to calculate the union of the re- 
gion segments, the random variables once defined can not be changed for different 
limit states. Since in the pressure profile the different pressure levels are presented, 
pressure should not be defined as a random variable. Therefore only two random 
variables are involved in gfun.dat - modulus of elastisity E and coefficient of thermal 
expansion a. The gfun.dat and various gfun.mov files Me shown in Appendix F. 

There are two ways of quantifying the effect of progressive damage on critical 
failure. The first is simply to compute the variation of critical failure probability 
with respect to progressive damage. This is shown in Fig. 26 for various pressure 
levels. 

An alternate way is to compute the progressive damage on overall critical failure 
probability. If each critical failure limit sate segment for each progressive damage 
gives the event F,, then the overall critical failure probability is P(U" =1 P \ where n 
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is the number of limit states ( n = 7 in this case ). 

The output files gfun.mov provide g-functions defining different limit-states for 
each level of damage as 

9i = Pi + £ (82) 

i'=i 

From Table 8, it is cleax that temperature variations do not have any significant 
effect on the structured failure probability. Therefore, only pressure variations are 
considered as follows: 

The parameters provided by gfun.mov for the structural limit states correspond- 
ing to fluid pressure profile are as follows: 


Pressure 

(3 

c*i 

<*2 

2700.0 psi 

6.825696 

0.999887 

-0.015057 

2900.0 psi 

5.721588 

0.999886 

-0.015088 

3000.0 psi 

4.619647 

0.999886 

-0.015071 

3200.0 psi 

2.423429 

0.999887 

-0.015038 

3300.0 psi 

1.328801 

0.999887 

-0.015021 

3400.0 psi 

0.236641 

0.999888 

-0.014959 

3420.0 psi 

0.018494 

0.999887 

-0.015003 


The above data provides parameters for 7 g-functions. Using the above data 
to calculate the union of the region defined by a group of g-functions [Ang and 
Tang, 1984], we obtain the probability of structural failure involving the progressive 
damage caused by fluid pressure. The probability is defined by a lower and upper 
bounds, which in this case are both 0.4926211. The data and output files are shown in 
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Appendix G. Those bounds provide the overall structural failure estimation when the 
system experiences various levels of progressive damage. The fluid pressure changes 
in a range from 2700.0 psi to 3420.0 psi. 

The above two methods provide practical tools for multi-disciplinary system re- 
liability estimation using NESSUS. With multiple impositions of one mode on the 
other mode, a close approximation to the failure domain can be constructed, and the 
critical failure probability can be obtained through the union of the failure region 
defined by the various limit-states. 
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CHAPTER V 


CONCLUSIONS AND RECOMMENDATIONS 
Conclusion 

This report has demonstrated the application of equivalence concepts to the re- 
liability analysis of multi-disciplinary systems using NESSUS. A thermal-structural- 
fluid system is used to illustrate the proposed methodology. The analogous model 
is a very powerful tool to analyze the one-dimensional steady state problem in heat 
transfer and fluid mechanics by converting those models into a structural model. 
Then the NESSUS probability analysis program can be implemented and the precise 
system reliability can be evaluated. Both traditional and progressive system fail- 
ure probability methods using NESSUS provide practical tools for multi-disciplinary 
system reliability analysis. 

Recommendations for future research 

This research project demonstrated how the NESSUS program could be applied 
for reliability analysis of engineering systems involving different disciplines, such as 
structure, heat transfer and fluid mechanics. The current models are based on the 
condition of one-dimensional, steady-state for both heat transfer and fluid mechanics. 
More complex systems could be treated in the similar way. However, the scope of 
application of this methodology is largely dependent on the ability of NESSUS/FEM 
to deal with problems in different disciplines under more complicated situations, for 
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instance, a thermal or a flow model in two or three dimensions and non-steady state 
conditions. For more sophisticated systems, the need for separate FEM program 
may be inevitable. Either a new source code should be developed or the existing 
commercial softwares could be integrated into the program. Several users of NESSUS 
have already integrated its FPI module to other FEM analysis program such as 
ANSYS and NASTRAN. Nevertheless, the use of equivalent concepts helps to obtain 
a quick estimate of multi-disciplinary system reliability through the use of NESSUS. 

The idea of progressive damage imposition to quantify nonlinear system reliability 
effects has previously been applied to structural mechanics problems [1]. This study 
extends this concept to multi-disciplinary systems. This appears to be a practiced 
methodology for system reliability analysis when failure modes (even in different 
disciplines) have physical relationship with each other. The methodology should be 
pursued further for application to other, more complicated engineering systems. 
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Appendix A 


FEM FILES FOR HEAT TRANSFER OF A HEAT EXCHANGER 


(Refer to Chapter III, Fig.16) 


NESSUS/FEM input file 


•xp2fam.dat 


Had Hay 14 20 1 37 *04 1997 


C . . . H1AT TRANSFER FOR A RICRICERANT PIPE 
C 

c 

c 

C ... PARAMETER DATA 
C 

20 
0 


•BOW 

•COM3 

•DISP 

•MODI 

•DOTL 

•SLIM 

91 

•PORC 
4 PTRT 
•PRIM 


•DTD 


MODEL DATA 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1 

3 

3 


0.0000 

10.0000 

10.0000 

20.0000 

20.0000 

30.0000 


0.0000 
0 . 0000 
0.0000 
0.0000 
0.0000 
0.0000 


0.0000 
0.0000 
0.0000 
0.0000 
0 . 0000 
0.0000 


0.0000 
0.0000 
0.0000 
0 . 0000 
0.0000 
0.0000 


1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 


( 

1 

0.0000 





I 

3 

0.0000 





1 

1 

0.0000 





1 

2 

70.0000 





1 

3 

0.0000 





1 

4 

0.0000 





1 

5 

0.0000 





1 

< 

0.0000 





•BEAM 

1 






1 

t 

1.0000 

1.0000 




•ITERA 

0 

4 





40 

0.0500 





•PROP 

1 

98 





2.51-4 

2 

1.0000 

69.43667 

0.3000 

0.0000 

3 

4 

1.0000 

79.9230 

0.3000 

0.0000 

2 . 51-4 

S 

6 

1.0000 

39955.3476 

0.3000 

0.0000 

2 . 3E-4 

•FORC 

< 

2 

-241.64 





•PRIN 

TOTA 

NODI 
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NESSUS/FEM output file 


axp2fam.out 

i 


W«d May 


14 20 ; 50 1 14 1997 


1 


NN NN EEEEEFEI SSSSSSS 

NNN NN El S3 

NN N NN EC 38 

NN H NN EEEZEt S3SSSS 
NH N NH EC SS 

NH NNN EC SS 

NH NN EEEEEEEI SSSSSSS 


SSSSSSS to uu 

S3 UU TO 

SS TO TO 

SSSS5S TO TO 

SS TO TO 

SS TO TO 

SSSSSSS TOTOTO 


SSSSSSS 

SS 

SS 

SSSSSS 

SS 

SS 

SSSSSSS 


NKNNN HNNH 

H H 

NN NNN NN 
NNNNNNN 
NNN NN NN 
N N N 

N N N 

W N N 


FFFFFFFF EEEEEEEC MM KM 
FF ES MMX KMM 
FF R KM M N MM 
FFFFFF CITE EE MM N MM 
FF R KM MM 
FF SC KM KM 
FF ESBEEEES KM KK 


DATS: 14- 3-1997 20:30 - LEVEL (.0 (109) - DATED Jun«. 30 , 1992 


C . . . HEAT TRANSFER FOR A RECRIOERANT PIPE 


•* PERTURBATION MO. 4 ** 


SUMMARY OF PERTURBATION DATA 


VARIABLE NO. 
PERTURBED VALUE 
VARIABLE NO. 
PERTURBED VALUE 

VARIABLE NO. 
PERTURBED VALUE 
VARIABLE NO. 
PERTURBED VALUE 


1 MEAN VALUE 

0.000000*00 

2 MEAN VALUE 

0.000000*00 

3 KEAN VALUE 

o.oooood*oo 

4 KEAN VALUE 

0.000000*00 


0.000000*00 STANDARD DEVIATION 
0.000000*00 STANDARD DEVIATION 
0.000000*00 STANDARD DEVIATION 
0.000000*00 STANDARD DEVIATION 


0 . 000000*00 
0.000000*00 
0.000000*00 
0.000000*00 


PERTURBATION ASSEMBLY, INCA. 0 PERT. 4 ITER. 0 CPTIKE- 0.31 SEC 


••• PERTURBATION 4 HAS NO EFFECT ON THIS ANALYSIS 


\ 

DATABASE UPDATE, INCA. 0 PERT. 4 ITER. 0 CPTIKE- 0.12 SEC 
1C . . . HEAT TRANSFER FOR A REGRIGERANT PIPE DATE 

i 14- 5-1997 20:50 RESULTS PAGE i 1 


TOTAL DISPLACEMENTS INCREMENT 0 TOTAL TRANSIENT TIKE 0.000000*0 

0 


NODE 
COMP. « 

COMP. 1 

COHP. 2 

COMP. 3 

COHP. 4 

COMP. 5 


1 

. 0OOO0D*00 

0 . 000000*00 

0 . 700000*02 

0.000000*00 

0.000000*00 

0.000000*00 

0 

2 

.000000*00 

0.000000*00 

0 . 1S296D*02 

0.000000*00 

0.000000*00 

0.000000*00 

0 

3 

.000000*00 

0.000000*00 

0 . 352980*02 

0 . 000000*00 

0.000000*00 

0.000000*00 

0 

4 

. 00000D*00 

0 . 00000D*00 

0 . 506 14D*01 

0.000000*00 

0.000000*00 

0.000000*00 

0 

5 

.000000*00 

0 . OOOOOD*00 

0 . 506 HD* 01 

0.000000*00 

0.000000*00 

0.000000*00 

0 

6 

. 00Q00D*00 

0.000000*00 

0.500090*01 

0.000000*00 

0.000000*00 

0.000000*00 

0 

1C . . . HEAT 
) DATE: 14- 

TRANSFER FOR A 
5-1997 20:50 

REGRIGERANT I 
PAGE: 11 

PIPE 


VERSION 4.0(109 


END OF INCREMENT, I VCR. 0 PERT. 0 ITER. 0 CPTIKE- 0.34 SEC 

1C ... HEAT TRANSFER FOR A REGRIGERANT PIPE VERSION 4.0(109 

> DATE: 14- 5-1997 20:50 PAGE: 12 


35 
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Appendix B 


FEM FILES FOR FLUID FLOW IN A DUCT 

(Refer to Chapter III, Fig. 18 ) 

NESSUS/FEM input file 


flowfaal.dat Sat Jun 28 22:28:43 1997 


1 


TEH 

C ... FLOW IN A CONSTANT DIAMETER PIPE (with A/L-0.1) 
C 
C 
c 

c - . . PARAMETER DATA 
C 

•BOON 20 

•CONS 0 

•orsp 

•NODE 2 

•EL EM 1 

38 

•FORC I 

•PERT 4 4 

• PR IN 

•BEAM 

•END 


C'"* MODEL DATA 


•COOR 

1 

2 

'EL. EM 
1 

'SOUN 

2 

2 

I 

1 

1 

1 

1 

1 

'BEAM 

1 

•ITEAA 

40 

• PROP 

1 

'FORC 

2 

•PRIM 


0.0000 

0.0000 

98 

1 2 

1 0.0000 

3 0.0000 

1 0.0000 

2 1138.7 

3 0.0000 

4 0.0000 

5 0.0000 

« 0.0000 

2 1.0000 

0 4 

0.0500 

98 

2 1.0000 

2 -240.25 


0.0000 

10.000 


1.0000 


14.8960 


0.0000 0.0000 

0.0000 0.0000 


0.3000 0.0000 


TOTA NODE 

'END 


0.0000 

0.0000 


2.SE-4 


1.0000 

1.0000 




72 


NESSUS/FEM output file 


flowfaml .out 


Sat Ju n 28 22:27:00 1997 


) CATS : 2®- 6-1997 22:2® PAGE : 9 


4 


•• PERTURBATION NO. 4 ** 


SUMMARY OF PERTURBATION DATA 


VARIABLE NO. 
PERTURBED VALUE 
VARIABLE NO. 
PKRTURBSD VALUE 
VARIABLE NO. 
PERTURBED VALUE 
VARIABLE NO. 
PERTURBED VALUE 


1 KEAN VALUE 
0. Q0000D*00 

2 MEAN VALUE 
0 . OGOOOD*O0 

3 MEAN VALUE 
0.000000*00 

4 KEAN VALUE 
0.000000*00 


O.OOOOOD*CO 

0.000000*00 

0.000000*00 

0.000000*00 


STANDARD DEVIATION 
STANDARD DEVIATION 
STANDARD DEVIATION 
STANDARD DEVIATION 


0.000000*00 
0.000000*00 
0 . OOOO0D*OO 
0.000000*00 


PERTURBATION ASSEMBLY, INCR. 0 PERT. 4 ITER. 0 CPTIME- 0.27 SEC 


■ ' * PERTURBATION 4 HAS NO EFFECT ON THIS ANALYSIS 


DATABASE UPDATE, INCR. 0 PERT. 4 ITER. 0 CPTIKE- 0.27 SEC 
1C . . - FLOW IN A CONSTANT DIAMETER PIPE (With A/L-O.U DAT 

E: 21- C-1997 22:25 RESULTS PACE: I 


00 


TOTAL DISPLACEMENTS 


INCREMENT 0 TOTAL TRANSIENT TIME 0.000000* 


NODE COMP. I COMP , 2 COMP. 3 

COMP. 5 

1 0.000000*00 0 . 1S3870* 04 0.000000*00 

0 . 0OQ0CD*O0 

2 0.000000*00 0 . 16774D*04 O.OO00OD»OO 

0.000000*30 

1C . . . FLOW IN A CONSTANT DIAMETER PIPE (with A/L-0.M 
9> DATS: 30- 4-1997 22:25 PACE; 10 


COMP. 4 


COMP. 5 


0.000000*00 0.000000*00 
0.000000*00 ^ 0.000000*00 

VERSION (.0(10 


END OF INCREMENT, INCR. O PERT. 0 ITER. 0 CPTIME- 0.29 SEC 

1C . . rLOW IN A CONSTANT DIAMETER PIPE (with A/L-0.11 VERSION (.0(10 

9) DATE: 21- (-1997 22:2( PACE: 11 


45 


STOP DUE TO END OF INPUT PILE 
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Appendix C 


PFEM FILES FOR STRUCTURAL RELIABILITY 

(Chapter IV, Data is shown in Fig. 16 of Chapter III) 

NESSUS/PFEM input and output files for single Z-level 


•xp3spfaa.dat Tu« Jun 24 10s 13: 00 1997 1 


• PEEK 

3 

-4.500 

0.95 

C . . „ PROBABILISTIC ANALYSIS POR STRUCTURAL RELIABILITY OR EEP2 

3 

-4.000 

0.95 

•MVDEPINE 

4 

-3.500 

0.95 

•COND 0 

3 

-3.000 

0.95 

•DATATYPE 0 

6 

-2.700 

0.95 

•RESPTYPB 3 

7 

-2.400 

0.95 

•COKP 3 

9 

-2.100 

0.93 

•NODI 44 

9 

-1.900 

0.95 

•PERT 3 

10 

-1.300 

0.95 

1.2,3 

11 

-1.200 

0.93 

•RANVAR 3 

12 

-0.900 

0.93 

1,2.3 

13 

-0.400 

0.95 

•DO 

14 

-0.300 

0.93 

•ZPDEPOT 

15 

0.000 

0.95 

•COMPUTATIONAL 1 3 

19 

0.300 

0.95 

1,2,3 

17 

0.400 

0.95 

*EHD 

19 

0.900 

0.95 

•RVDEFIN1 

19 

1.200 

0.95 

•DEPOT 1 

20 

1.300 

0.95 

PI 

21 

1.900 

0.95 

1377.4 313.41 NORMAL 

22 

2.100 

0.95 

PRESSURE 

23 

2.400 

0.95 

30 31 1.0 

34 

2.700 

0.95 

•DEPOT 2 

25 

3.000 

0.95 

EMOD 

24 

3.500 

0.93 

0 . 1700000B*09 0 . 17000008*07 NORMAL 

27 

4 .O0> 

0.95 

PROP 153 

29 

4.300 

0.95 

1 59 0.0 1.0 0.0 0.0 0.0 

29 

5.000 

0.95 

•DEPINE 3 

30 

-5.000 

0.445 

ALPA 

31 

-4.500 

0.445 

0. 9300000E-05 0 . 1900000K-05 NORMAL 

32 

-4.000 

0.443 

PROP 153 

33 

-3.500 

0.445 

1 31 0.0 0.0 0.0 1.0 0.0 

34 

-3.000 

0.445 

•PERT 1 

33 

-2.700 

0.445 

1 0.1 

34 

-2.400 

0.445 

•PERT 2 

37 

-2.100 

0.445 

2 0.1 

39 

-1.000 

0.4(5 

•PERT 3 

39 

-1.300 

0.445 

3 0.1 

40 

-1.200 

0.445 

•END 

41 

-0.900 

0.445 

•END 

42 

-0.400 

0.445 


43 

-0.300 

0. 445 

•n k 

44 

0.000 

0.445 

C . . . REFRIGERANT PIPE UNDER INSIDE PRESUME 

45 

0.300 

0.445 

C 

44 

0.(00 

0.(45 

c 

47 

0.900 

0.445 

c 

49 

1.200 

0.445 

C . . . PARAMETER DATA 

49 

1.500 

0.(45 

C 

30 

1.900 

0.445 

•CONS 0 

31 

2.100 

0.445 

•DISP 

52 

2.400 

0.445 

•NODE 39 

33 

2.700 

0.445 

•ELEN 29 

54 

3.000 

0.445 

153 

35 

3.500 

0.445 

•BOUN 20 

34 

4.000 

0.645 

•PRESSURE 

57 

4.500 

0.665 

•TEMPERATURE 

59 

5.000 

0.665 

•LOUB 2 1.3 1 

•ELEN 

153 


•PERT 10 10 

1 

1 30 

31 2 

PRESSURE 

2 

2 31 

32 3 

PROPERTIY 

3 

3 32 

33 4 

PROPERTY 

4 

4 33 

34 5 

•PRIN 

3 

5 34 

35 ( 

•END 

4 

4 35 

36 7 

C*«»* MODEL DATA 

7 

7 34 

37 9 

•COOR 

9 

9 37 

31 9 

1 -3.000 0.95 

* 

9 39 

39 10 
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•xp2sp£< 


2 


«a.d*t Tu« Jun 24 18tl3t00 1997 


10 

10 

39 

40 

11 

11 

11 

40 

41 

12 

12 

12 

41 

42 

13 

13 

13 

42 

43 

14 

14 

14 

43 

44 

15 

13 

IS 

44 

43 

14 

1( 

1C 

43 

4C 

17 

17 

17 

44 

47 

IS 

11 

11 

47 

41 

19 

19 

19 

41 

49 

30 

20 

20 

49 

50 

21 

21 

21 

50 

51 

22 

22 

22 

51 

52 

23 

23 

23 

52 

53 

24 

24 

24 

33 

34 

23 

23 

25 

34 

35 

24 

2S 

2C 

33 

SC 

37 

27 

27 

5C 

57 

31 

21 

2B 

57 

51 

39 

•PRES SUM 




30 

SI 

1577. 

4 


•TOT? ERA TOM 




1 

29 

4.06 



30 

31 

4.00 



•BOUH 





IS 

1 

0.0 



44 

1 

0.0 



•ITERA 

1 3 




40 

0.0300 





1 51 0.1000t*01 0 . 17008*08 0.340 0.9500E-0S 0.3UI «.0t*0« 

•PHD* 

STM MODI 



HEAT EXCKACER OP EXP 2 WITH DATASETS 
•RVNUM 3 
•orcwcTiaK i 
•DATASETS 4 
•METHOD 1 
•PRINTOPT 0 
•ANALTYPE 1 
•END 

•ZLZVSLS 1 
0. 1000001*04 
•END 


DESIGN SENSITIVITIES 

TAYX-OR series expansion or THE form 
G ■ A0 ♦ A1 *X1 ♦ A2*X2 ♦ ... ♦ AN*XN 


•xp2*p£«A. out Tu« Jun 24 18tl3t42 1997 3 

A0 IS THE CONSTANT TERM 

Al. A3, ... , AN AM THE DESIGN SENSITIVITIES 
XI, 12 XN AM THE RANDOM VARIABLES 


TAYLOR SERIES COEFTICIENTS 
A COEPTS . VALUE 


RANDOM EXPANSION NORM. DESIGN 

VARIABLE POINT (KPP) SENSITIVITY 


0 0. 2098461*04 

1 0 . 2333341*01 

2 -0.1233211-03 

3 0.1417301*07 


1 0.1377401*01 0.(333431*00 

2 0.1700001*01 -0. 341BB9E*00 

3 0. 9300001-05 0.27(7301-02 


PERTORMING PROBABILISTIC ANALYSIS WITH PPI 


Z 

0 . 800000001*04 


CDF RESULTS 

U PROBABILITY 

0 . 3(2237(3E*01 0 . 99999999E*00 
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NESSUS/PFEM input and output file for CDF 


•xp2scdf.dat Tu» Jun 24 19il9t37 1997 1 


•Fra 

2 

-4.300 

0.95 

C ... PROBABILISTIC ANALYSIS FOR STRUCTURAL RELIABILITY OP KP2 

3 

-4.000 

0.93 

•MVDEFINE 

4 

-3.300 

0.93 

•CCND a 

3 

-3.000 

0.95 

•DATATYPE 0 

4 

-2.700 

0.95 

•RXSPTYPI 3 

7 

-2.400 

0.93 

•COMP 3 

• 

-2.100 

0.95 

•NODI 44 

5 

-1.800 

0.95 

•POT 3 

to 

-1.300 

0.93 

1.2,3 

11 

-1.200 

0.93 

•RANVAR 3 

12 

-0.900 

0.93 

1.3,3 

13 

-0.400 

0.93 

'DID 

14 

-0.300 

0.93 

•IFDEFIN* 

13 

0.000 

0.93 

•COMPUTATIONAL l 3 

16 

0.300 

0.93 

1.2.3 

17 

0.400 

0.93 

•END 

18 

0.900 

0.95 

•RVDEFWE 

19 

1.200 

0.93 

•define i 

20 

1.300 

0.93 

PI 

21 

1.800 

0.93 

0 . 137740001*04 0.31548001*03 NORMAL 

22 

2.100 

0.93 

PRESSURE 

23 

2.400 

0.93 

30 31 1.0 

24 

2.700 

0.93 

•DEFINE 2 

23 

3.000 

0.93 

EFfOD 

24 

3.300 

0.93 

0 . 17000001401 0.1700000**07 NORMAL 

27 

4.00<F 

0.95 

PROP 133 

28 

4.300 

0.95 

1 SB 0.0 1.0 0.0 0.0 0.0 

29 

5.000 

0.95 

•ORPINE 3 

30 

-5.000 

0.445 

ALFA 

31 

-4.500 

0.445 

0 . 93000001-05 0.19000001-03 NORMAL 

32 

-4.000 

0.445 

PROP 133 

33 

-3.500 

0.445 

1 38 0.0 0.0 0.0 1.0 0.0 

34 

♦3.000 

0.443 

•PERT 1 

35 

-2.700 

0.445 

1 0.1 

34 

-2.400 

0.445 

• PERT 2 

37 

-2.100 

0.445 

2 0.1 

38 

-l ♦ BOO 

0.445 

•PERT 3 

39 

-1.500 

0.443 

3 0.1 

40 

-1.200 

0.443 

•END 

41 

-0.900 

0.445 

•DID 

42 

-0.400 

0.445 


43 

-0 . 300 

0 . 445 

•FEX 

44 

0.000 

0.445 

C ...REFRIGERANT PIPE ONDER INSIDE PRESURE 

43 

0.300 

0.443 

C 

44 

0.400 

0.443 

C 

47 

0.900 

0.443 

C 

48 

1.200 

0.443 

C . . . PARAMETER DATA 

43 

1.500 

0.443 

C 

50 

1 . 800 

0.645 

•CONS 0 

31 

2.100 

0.665 

•DISP 

52 

2.400 

0.463 

•NODE 38 

53 

2.700 

0.465 

•ELXM 28 

54 

3.000 

0.663 

153 

55 

3.500 

0.663 

•BOUN 20 

54 

4.000 

0.663 

•PRESSURE 

57 

4.500 

0.663 

•TEMPERATURE 

58 

5.000 

0.665 

•LOOT 2 1 ' 3 1 

•ILEX 

153 


•PERT 10 10 

1 

1 30 

31 2 

PRESSURE 

2 

2 31 

32 3 

PROPERTIY 

3 

3 32 

33 4 

PROPERTY 

4 

4 33 

34 3 

•PRIM 

3 

5 34 

35 6 

•END 

4 

4 35 

36 7 

C ,f ** MODEL DATA 

7 

7 36 

37 I 

•COOR 

8 

8 37 

38 9 

1 ' ' -3.000 ~ 0.53 

9 

9 38 

39 10 

- r ’j-. -* 
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•xp2aedf.dat 


Tu* Jua 24 19 : 19 i 57 1997 


2 


10 

10 

39 

40 

11 

11 

11 

40 

41 

12 

13 

13 

41 

43 

13 

13 

13 

42 

43 

14 

14 

14 

43 

44 

15 

13 

IS 

44 

43 

16 

16 

If 

43 

46 

17 

17 

17 

46 

47 

IB 

18 

18 

47 

48 

19 

19 

19 

41 

49 

20 

30 

30 

49 

50 

21 

31 

31 

30 

51 

32 

33 

33 

31 

53 

33 

33 

33 

S3 

S3 

34 

34 

34 

33 

34 

25 

35 

35 

34 

35 

26 

36 

36 

35 

36 

27 

37 

27 

36 

37 

28 

31 

28 

37 

38 

29 


•PRESSURE 

30 51 0.15774001*04 

•TEMPERATURE 

l a> i.o< 

30 5« 4.00 

•»OUN 

13 1 0.0 

44 1 0.0 

•ITIM 1 3 

40 0.0500 

•PROP 133 

1 51 0.10001*01 0.17008*08 0.340 0.95001-05 0.3168 6.08*06 

•PRIM 

•TR* IKJD* 

•KKD 



•PPI 

HEAT ECCHAGER OP IXP2 WITH DATASETS 
•RVNUW 3 
•UFUNCTION 1 
•DATASETS 4 
•METHOD 1 
•FRINTOPT 0 
•ANALTYPI 0 
•END 
•END 


DESIGN SENSITIVITIES 

TAYLOR SERIES EXPANSION OP THE FORM 

O ■ A0 * A1*X1 ♦ A3*X2 * ... + AN*XN 
WHERE i 


•xp2«cdf . out Tu# JUn 34 18t38:49 1997 3 

A0 IS THE CONSTANT TERN 

Al, A3, ... , AN ARE THE DESIGN SENSITIVITIES 
XI, X2 , ... , XN ARE THE RANDOM VARIABLES 


TAYLOR SERIES COEFFICIENTS 
A COEPPS. VALUE 


RANDOM EXPANSION NORM. DESIGN 

VARIABLE POINT <KPP> SENSITIVITY 


0 0.309<4<B*04 

1 0.3333341*01 1 
3 -0.1233311-03 3 
3 0.1667501*07 3 


0. 1377408*04 0.6353438*00 
0.1700008*06 -0.3616696*00 
0 . 9 500008-03 0.376730E-03 


PERFORMING PROBABILISTIC ANALYSIS WITH PPI 


-0 . 306899416*03 
0. 551499691*03 
0.133389931*04 
0.20*939851*04 
0.2854697*1*04 
0.431549641*04 
0.513089571+04 
0 . 59 1629 50E+04 
0. 668169438+04 
0 . 744709 J5E+04 


CDF RESULTS 
U 

-0 . 5 10000006*01 
-0.41000000E+01 
-0.310000006*01 
-0 . 210000008*01 
-0.110000001*01 
0 . 900000008*00 
0.190000008*01 
0 . 390000006*01 
0 . 39000000E*0 I 
0 . 49000000E+01 


PROBABILITY 
0 . 170123 3 IB-06 
0.306687168-04 
0 . 96767132E-03 
0. 178 643 57E-01 
0.135666108*00 
0.81593991E+00 
0.971283 51E+00 
0. 998134128+00 
0. 999951888*00 
0 . 99999952 E+00 


ITER. NO. 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
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Appendix D 


PFEM FILES FOR THE HEAT TRANSFER RELIABILITY 


(Chapter IV, data is shown in Fig. 16 in Chapter III) 

NESSUS/PFEM input and output files for single Z-level 


•xp2tompftmx.dat Sun Jun 22 01:05i39 1997 1 

•pfd* 

C . . . PROBABILISTIC ANALYSIS MU TEMPERATURE FAILURE OF EXP2 I -LEVEL 
•HVDEFINI 

•COND 0 

•DATATYPE 0 
•RXSPTYPE I 
•COK? 2 

•NOD* < 

•FEW 3 

1.3.3 

•RANVAR 3 

1.3.3 

•END 

•ZFDEFIN1 

•COMPUTATIONAL 1 3 


•END 

•KVDEFXNE 

•DEFINE 1 

EM0D3 

0 .39935(7(0001405 0 . 3S95$47«O00E*04 NORMAL 

PROF SI 

5 I 0.0 1.0 0.0 0.0 0.0 

•DEFINE 2 

EHOD3 

0 . 5117310008*02 0 . 58I758000E*01 NORMAL 

PROF SI 

3 I 0.0 1.0 0.0 0.0 0.0 

•DEFINE 3 

EMODI 

0 .7853511*03 0.7I53SIE+01 NORMAL 

FROF SI 

1 3 0.0 1.0 0.0 0.0 0.0 

•PERT 1 

1 0.1 

•PERT 2 

2 0.1 

•PIET 3 

3 0.1 

•END 

•END 

C < 

•FEK 

C . . . HEAT TRANSFER FOR A REQAIGERANT PIPE EXP2 

C 

C 

c 

C . . . PARAMETER DATA 
C 

•BOON 20 

•CONS 0 

•DISP 

•NODE I 

•DUEL 4 

•ELEN 3 

SI 

•FORC 1 

•PERT 10 10 / 

PROPERTIES 

PROPERTIES 

PROPERTIES 

•FRIN 

•BEAM 

•END 

C'*'» MODEL DATA 

•COOR 

1 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 
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2 

3 

« 

5 

C 

•ELEX 

1 

2 

3 

*DUPL 

3 

5 

•BOUN 

€ 

6 
1 
1 
1 
1 
1 
1 

‘BEAM 

1 

•ITERA 

40 

•PROP 

l 

3 

5 

*FORC 

< 

•PR IK 




0.0000 10.0000 
0.0000 10.0000 
0.0000 20.0000 
0.0000 20.0000 
0.0000 30.0000 

1 2 

3 4 

J < 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


0.0000 

0.0000 

0.0000 

70.0000 

0.0000 

0.0000 

0.0000 

0.0000 


1.0000 1.0000 

3 

L0500 

8 > 

I 1.0000 78.3390 

l 1.0000 58.1751 

S 1.0000 39953 . 475 


0.3000 

0.3000 

0.3000 


2 


-241.64 


TOTA MOD* 


0.0000 
0.0000 
0.0000 
0 . 0000 
0.0000 


0.0000 
0.0000 
0 . 0000 


0.0000 

1.0000 

0.0000 

1.0000 

0.0000 

1.0000 

0.0000 

1 . 0000 

0.0000 

1.0000 


2 . 51-4 
2 . 51-4 
2. St-4 


•DID 

***** 

*FPI 

REFRIGERANT PIP* RELIABILITY WITH Z-LEVELJT < 5*P1 
•RVNUM 3 
•GFUNCTION 1 
•DATASETS 4 
•METHOD 1 
♦PRINTOPT 0 
•ANALTYPE 1 
•DID 

•ZLEVELS 1 

3.00000 

•END 


•xp2t«jnpf tmz . out Sun Jun 22 01x06:16 1997 3 

A0 IS THE CONSTANT TERM 

Al, A3, ... , AM ARE TK* DESIGN SENSITIVITIES 
XI, X2 XN ARE THE RANDOM VARIABLES 


TAYLOR SERIES COEFFICIENTS RANDOM EXPANSION WORM. DESICH 

A COEFFS. VALUE VARIABLE POINT (HPPJ SENSITIVITY 


0 -0 . 7375041*02 

1 0.151374E-03 

2 0 .4971581+00 

3 0.3917651*00 


1 0.3995551*05 0.841489E-03 

2 0.5887588*02 0.5710«8B*00 

3 0.7833981*02 0.4210901*00 


PERFORMING PROBABILISTIC ANALYSIS WITH FPI 


Z 

0.300000001*01 


CDF RESULTS 
U 

0.134037131*01 


PROBABILITY 
0 . 909 921 40E«00 


ITER. NO. 
0 
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NESSUS/PFEM input and output files for CDF 


•xp 2twnpfam.dat Thu Apr 6 17tl4 *48 1993 

•PT1M 

C . . . PROBABILISTIC CCDP) ANALYSIS FOR TEXPERATORE FAILURE OF UM 
•KVDETIHE 

•COKD 0 

•DATATYPE 0 

•Msmri 1 
•cow a 
•wool i 
•tot a 

1 . 2.3 

•RAMVAR 3 

1.3.3 

•BO 

•ZFDSPZXB 

•COWOTATIONAL 1 3 

1.2.3 


•RVDEFXNX 

•DEFINE 1 

EM0D3 

0 . 39953 474000E* 05 0.399554790001*05 WORKAL 
PROP 91 

3 « 0.0 1.0 0.0 0.0 0.0 

•DEPTHS 2 

D40D2 

0.3917590001*02 0.5997310001*01 NORMAL 
PROP 99 

3 4 0.0 1.0 0.0 0.0 0.0 

•DSP INI 3 

ZMOD1 

0.7953991*02 0.7933991*01 NORMAL 
PROP 99 

1 a o.o i.o o.o o.o o.o 

•PSRT 1 

1 0.1 

•PSRT 2 

2 0.1 

•PSRT 3 

3 0.1 

• BID 

• 1ND 

C * ...♦•♦.it.* * 

*rrM 

C . . . HEAT TRANSFER. FOR A REGRIGERANT PIPS EXP 2 
C 

c 

c 

C . . . PARAMETER DATA 
C 

•10UN 20 

•CONS 0 

•DISP 

•NODI C 

•DCJPL 4 

•SUM 3 

99 

•FORC 1 

•PSRT 10 10 ' 

PROPERTIES 
PROPERTIES 
PROPERTIES 
•PR IN 
•BEAM 
•END 

C* ,B * MODEL DATA 
•COOR 

1 0.0000 0.0000 0.0000 0.0000 0.0000 1 


1 


0000 
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2 

0.0000 

10.0000 

0.0000 

0.0000 

0.0000 

1 . oooo 

3 

0.0000 

10.0000 

0.0000 

0.0000 

0.0000 

1.0000 

4 

0.0000 

20.0000 

0.0000 

0.0000 

0.0000 

1.0000 

5 

0.0000 

20.0000 

0.0000 

0.0000 

0.0000 

1.0000 

6 

0.0000 

30.0000 

0.0000 

0.0000 

0.0000 

1.0000 


•ELEM 

1 

2 

3 

*DUPL 

3 

5 

*BOUN 


0.0000 
0.0000 
0 . 0000 
70.0000 
0.0000 
0.0000 
0.0000 
0.0000 


•ITERA 

40 

•PROP 


3 

•FORC 

6 

•PR IN 


•END 

C 

•FPI 


0000 


1.0000 

1.0000 

1.0000 


TOT* NODE 


78.5398 

58.8758 

35855.478 


0.3000 

0.3000 

0.3000 


0.0000 

0.0000 

0.0000 


2.5E-4 

2.5B-4 

2.3B-4 


HEAT EXCHAGER OF EXP2 WITH DATASETS 
•RVNUX 3 
•GFUNCTION 1 
•DATASETS 4 
•METHOD 1 
•PRINTOUT 0 
•ANAL TYPE 0 
•END 
•END 


•xp2tarapfam.out Mon May 19 22:40:36 1997 3 

A0 IS THE CONSTANT TERM 

Al, A2 AN ARE THE DESIGN SENSITIVITIES 

XI, X2 XN ARE THE RANDOM VARIABLES 


TAYLOR SERIES COEFFICIENTS 
A COEFFS. VALUE 


RANDOM EXPANSION NORM. DESIGN 

VARIABLE POINT (MPP1 SENSITIVITY 


0 -0 . 737506E+02 

1 0. 151374E-05 1 

2 0 . 697 158E + 00 2 

3 0 . 39 1765E+0Q 3 


0 . 399555E+03 0.8414B9E-03 
0 . 588758E+02 0 . 571068E»00 
0. 7B5J98E+02 0 . 42 8090E*00 


PERFORMING PROBABILISTIC ANALYSIS WITH FPI 


Z 

-0 . 28039092E+02 
>0 . 2290894 1B+02 
-0. 17778790E»02 
-0. 1264B639E+02 
-O.7518408OE*O1 
0.2741 8 13 8E*01 
0 . 787196485*01 
0 . 1300211 6E*02 
0 . 10132267B+O2 
O.23262410B+O2 


CDF RESULTS 
U 

-0 . 51000000E+01 
-0 . 4100000QE+Q 1 
-0. 31000000E+01 
-0 . 21000000E+01 
-0 . 11000000E+01 
0 . 90000000E*00 
0 . 190000008*01 
0 . 29 0000008*01 
0 . 39000000E* 01 
0 . 490000008*01 


PROBABILITY 
0. 17012231E-06 
0 . 206607 16E-04 
0.967 6712 2E-03 
0 . 178643 57E-01 
0 . 1 356661 OE*00 
0 . 81593 991E*00 
0 . 971 283 51E*00 
0.998134128*00 
0 . 99995 1B8E*00 
0 . 999999528*00 


ITER. NO. 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
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Appendix E 


PFEM FILES FOR RELIABILITY OF FLUID FLOW 


(Chapter IV, data is show in Fig.18 of Chapter III) 

NESSUS/PFEM input and output files for a single Z-level 


flowrpftaz.dat Sat Jun 28 22:36:52 1997 1 

•prtx 

C PROBABILISTIC ANALYSIS FOR FLOW FAILURE 

•XVDEFINE 

•COND 0 

* DATATYPE 0 
•RESPTYPE 1 
•COM? 2 

•NODE 2 

•PERT 2 

1.3 

•RANVAR 2 

1.2 

•DO 

•ZFLEFINK 

•COMPUTATIONAL l 2 


•END 

•RVDEFINI 

♦DEFINE 1 

EMOD 

14.196 1.4*96 NORMAL 

PROP 9* 

1 2 0.0 1.0 0.0 0.0 0.0 

•DEFINE 2 

VELO 

-350.0 35.0 NORMAL 

F0RC 

2 2 1.00 

•PERT l 
1 0.1 
•PERT 2 
2 0.1 


*UfD 

•END 

C •« 


•n> 

C ... FLOW IN A CONSTANT DIAMETER PIPE With A/L-0.U 

C 

C 

c 

C . . . PARAMETER DATA 
C 

•BOUN 20 

•CONS 0 

■CIS? 

•NODE 2 

•ELIM 1 

91 

•FORC 1 

•PERT 10 10 

VELOCITY 

FORCE 

•PAIN 

•BEAM 

•END 

c*»«* model data 


•COOR 

1 

2 

•ELEM 

1 

•BOON 


0.0000 0.0000 

0.0000 10.000 

99 

l 2 

1 0.0000 

3 0.0000 

1 0.0000 

2 1838.7 

3 0.0000 


0.0000 

0.0000 


0.0000 

0.0000 


0.0000 

0.0000 


X. 


1.0000 

1.0000 
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I 4 0.0000 

1 5 0.0000 

1 6 0.0000 

* BEAM 1 

1 2 1.0000 

•ITERA 1 2 

40 0.0500 

-PROP 99 

1 2 1.0000 

•FORC 

2 2 -350.0 

-PRDI 

TOTA MODE 


1.0000 


14.1960 


0.3000 


0.0000 


2 . 5S-4 




TPI 

PLOW FAILURE PROBABILITY WITH Z-LEVEL 1 P<1677 . 4pa i ) 
•RVNUM 2 
•GFUWCTION 1 

•DATASETS 3 

•METHOD l 

•PRINTOPT 0 
•ANALYTYPE 1 

•SOTO 

'ZLEVXLS 1 

1677.4 

•END 


DESIGN SENSITIVITIES 

TAYLOR SERIES EXPANSION OP THE FORM 

O • A0 * A1*X1 ♦ A2 *X3 * ... * AN*XN 

WHERE: 

A0 IS THE CONSTANT TERN 

A l, X2 AN ARE THE DESIGN SENSITIVITIES 

Xl \ X2 XN ARE THE RANDOM VARIABLES 


TAYLOR SERIES COEFFICIENTS 
A COEFFS. VALUE 


0 0 . 1603748+04 


RANDOM EXPANSION NORM. DESIGN 

VARIABLE POINT IHPPl SENSITIVITY 


Clowpf 6BX .out 


Sat Jun 28 22i37tl3 1997 3 


1 0.1577358*02 

3 0.6713311*00 


1 0. 148960E*03 0.5000008*00 

3 -0. 150000E+03 -0 . 500000E+00 


PERFORMING PROBABILISTIC ANALYSIS WITH PPI 


CDF RESULTS 

j U PROBABILITY 

0.167740001*04 0 . 33 174074E*01 0. 98670241E*00 
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NESSUS/PFEM input and outpout files for CDF 


flowcdf.dat Sat Jun 28 22:45:21 1997 1 

•PPTX 

C PROBABILISTIC ANALYSIS FOR FLOW FAILURE 

•NVDEFINI 

•CONS 0 

•DATATYPE 0 
•RZSPTYPE 1 
•COKP 2 

•NODE 2 

•PERT 2 

1 , 2 

•RANVAR 2 
1.2 
•END 
•ZFDEFCNK 

•COKPUTATrONAL 1 2 

1.2 

•end 

•RVDEFINX 

•DEFINE 1 

IMOD 

14.896 1.4896 NORMAL 

PROP 98 

1 2 0.0 1.0 0.0 0.0 0.0 

•DEFINE 2 

VELO 

-350.0 35.0 NORMAL 

FORC 

2 2 1.00 
•PERT l 
1 0.1 
•PERT 2 
2 0.1 
•DID 
•END 

C •*** 

•HR 

C ... FLOW rN A CONSTANT DIAMETER PIPE {with A/L-0.1) 

C 

C 

c L 

C . . PARAMETER DATA 
C 

■»wn 20 

•CONS 0 

•DISP 

•NODE 2 

•ELEN 1 

93 

•porc l 

•PERT 10 10 

VELOCITY 
FORCE 
•PR IN 
•BEAM 
•END 

C*»»* MODEL DATA 
•COOR 

1 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 

2 0.0000 10.000 0.0000 0.0000 0.0000 1.0000 

•ELAM 99 

11 2 
•BOUN 

2 1 0.0000 

2 3 0.0000 

1 1 0.0000 

1 2 1833.7 

1 3 0.0000 
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1 4 0.0000 

l 5 0.0000 

1 6 0.0000 

‘BEAK 1 

1 2 i .0000 1.0000 

•ITERA 1 2 

40 0.0500 

•PROP 98 

1 2 1.0000 14.896 0.3000 0.0000 2.5E-4 

•FORC 

2 1 -350.0 

•PRXN 

TOTA NOOK 

* * 

*FPI 

FLOW FAILURE PROBABILITY WITH A SERIES Z-LEVEL5 
•RVNUW 2 
•□FUNCTION 1 
•DATASETS 3 

•METHOD 1 

•PRINTOPT 0 
• AKALYTYPE 0 

•END 
•END 


DESIGN SENSITIVITIES 


TAYLOR SERIES EXPANSION OP THE FORM 

C « A0 ♦ A1 *Xl ♦ A2 * X2 * ... + AN *XN 
WHERE: 

A0 IS THE CONSTANT TERM 

Al, A2 AN ARE THE DESIGN SENSITIVITIES 

XI, X2 XN ARE THE RANDOM VARIABLES 


TAYLOR SERIES COEFFICIENTS RANDOM EXPANSION MOM4. DESIQ* 

A COKFFS. VALUE VARIABLE POINT (MPP) SlNSTTIVITY 


0 0.1603748*04 


flowcdf.out Sat Jun 28 22:46:36 1987 3 

1 0 . 15773 58*02 1 0.1489608*02 0.5000008*00 

2 0.671321**00 2 -0 . 350000E*03 -0.5000008*00 


PERFORMING PROBABILISTIC ANALYSIS WITH FPI 


CTr RESULTS 

z 0 PROBABILITY ITER. NO. 

0.143431668*04 -0 . 51000000E*01 0. 170122318-06 0 

0 . 14 575 3 658*04 -0 . 41000000E*01 0 . 206687 1 6E-04 0 

0 . 150075648*04 -0 . J1000000t*01 0 . 967671228-03 0 

0. 153397638*04 -0.210000008*01 0 . 17 8643578-01 0 

0.1567lJ61g*04 -0 . 110000008*01 0.135666106*00 0 

0. 16336359E*04 0 . 90000000E*00 0 . 81 593991E*00 0 

0.16668553**04 0 . 19000000E*01 0 . »7128351E*00 0 

0.170007568*04 0 . 2 9000000E*0 1 0.998134128*00 0 

0.173329558*04 0.390000008*01 0.999951838*00 0 

0.176651548*04 0.490000008*01 0.999999528*00 0 
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Appendix F 


GFUN.DAT AND GFUN.MOV FILES OF IMPOSITIONS 


(Chapter IV, data is show in Fig.16 of Chapter III) 

NESSUS/PFEM input gfun.dat and output gfun.mov files 


gfun.dat Tu« JUl 1 15:39:34 1997 1 


* PFEM 



10 

-1.500 

0.95 

C - . . PROBABILISTIC ANALYSIS FOR STRUCTURAL RELIABILITY 

11 

-1.200 

0.95 

•HVDEFINE 



12 

-0.900 

0.95 

*COND 

0 


13 

-0.600 

0.95 

•DATATYPE 

0 


14 

-0.300 

0.95 

•RXSPTYPE 

3 


15 

0.000 

0.95 

•COMP 

3 


16 

0.300 

0. 95 

•NODE 

41 


17 

0.600 

0.95 

•PERT 

2 


18 

0.900 

0.95 




19 

1.200 

0.95 

•RANVAR 

2 


20 

1.500 

0.95 

1,2 



21 

1. BOO 

0.95 

•END 



22 

2.100 

0.95 

•ZFDEFINE 



23 

2.400 

0.95 

•COMPUTATIONAL 

1 2 

24 

2.700 

0.95 

1,2 



25 

3.000 

0.95 

•END 



26 

3.500 

0.95 

•RVDEFINE 



27 

4.000 

0.95 

•DEFINE 

1 


28 

4.500 

0.95 

EMOD 



29 

5.000 

0.95 

0 . 1700000E+08 

0 . 17000008*07 NORMAL 

30 

-5-000 

0.665 

PROP 153 



31 

-4.500 

0.665 

1 58 0. 

0 1 . 

0 0.0 0.0 0.0 

32 

-4.000 

0.665 

•DEFINE 

2 


33 

-3.500 

0 665 

ALFA 



34 

-3.000 

0.665 

0.9500000B-05 

0. 1900000E-05 NORMAL 

35 

-2.700 

0.665 

PROP 153 



36 

-3.400 

0.665 

1 58 0.0 

0.0 

0.0 1.0 0.0 

37 

-3.100 

0.665 

•PERT l 



38 

-1.800 

0.665 

1 0.1 



39 

-1.500 

0.665 

•PERT 2 



40 

-1.200 

0.665 

2 0.1 



41 

-0.900 

0.665 

•END 



42 

-0.600 

0.665 

•END 



43 

-0.300 

0.665 

•FEH 



45 

0.300 

0.665 

C ...REFRIGERANT PIPE 

UNDER INSIDE FRESUR1 

46 

0.600 

0.665 

C 



47 

0.900 

0.665 

C 



48 

1.200 

0.665 

c 



49 

1.500 

0.665 

C . . . PARAMETER DATA 


50 

1.800 

0.665 

C 



51 

2.100 

0.665 

•CONS 0 



52 

2.400 

0.665 

•DISP 



53 

2.700 

0.665 

•NODE 58 



54 

3.000 

0.665 

■ELEM 28 



55 

3.500 

0.665 

153 



54 

4.000 

0.665 

•BOON 20 



57 

4,500 

0.665 

•PRESSURE 



58 

5.000 

0.665 

•TEMPERATURE 



•ELEM 

153 


•LOUB 2 

1 3 

1 

1 

1 30 

31 2 

•PERT 10 

10 


2 

3 31 

32 3 

PROPERTY 



3 

3 33 

33 4 

PROPERTY 



4 

4 33 

34 5 

•PR IN 



5 

5 34 

35 6 

•END 



6 

6 35 

36 7 

C**" MODEL DATA 


7 

7 36 

37 8 

•COOR 



8 

8 37 

38 9 

1 

-5.000 

0.95 

9 

9 38 

39 10 

2 

-4.500 

0.95 

10 

10 39 

40 11 

3 

-4.000 

0.95 

11 

11 40 

41 12 

4 

-3.500 

0 . 95 

12 

12 41 

42 13 

5 

-3.000 

0.95 

13 

11 42 

43 14 

5 

-3.700 

0.95 

14 

14 43 

44 15 

7 

-2.400 

0.95 

15 

15 44 

45 16 

8 

-2.100 

0.95 

14 

16 45 

46 17 

5 

-1.800 

0.95 

17 

17 46 

47 18 
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Ofun.dat Tu* Jul 1 15 1 39 1 34 1997 3 


11 

18 

47 

48 

19 

19 

19 

48 

49 

20 

20 

20 

49 

30 

21 

21 

21 

30 

31 

22 

22 

22 

51 

52 

23 

23 

23 

32 

53 

24 

24 

24 

53 

54 

25 

25 

25 

54 

35 

24 

24 

26 

33 

54 

27 

27 

27 

34 

57 

28 

29 

•PRESSURE 

28 

37 

59 

29 

30 

59 

3300. 

0 



•twpiraturi 

1 29 4 . 0< 

20 59 4.0 

•BOON 

IS 1 0.0 

44 l 0.0 

•ITIRA 1 2 

40 0.0500 

•PROP 153 

1 39 0. 10001*01 0 . 17001*09 0.340 0.95001-05 0.3149 4.01*04 

•miN 

STR1 NODI 

•DID 

C 

•FPI 

HEAT KXCKAGER OP EXP 2 WITH OATASETS 
•RVNIIM 2 
•OTONCTION 1 
•DATASETS 3 
•METHOD 1 
■PAINTOPT 0 
•ANALTYP1 1 
•END 

•ZLEVELS 1 

0.9000001*04 

•END 


gfun.mov 

Tu« 

Jul 1 15 1 28 1 37 1997 

1 

CONDITION 1 

NODI 

COMPONENT 

LAYER 


INCA 0 

44 

3 

1 



Z MEDIAN, MEAN, STD. DEV. 

0.4549371*04 0.4549371*04 0.2125241*03 

NUMBER OP PROBABILITY LEVELS 1 

LVL, Z0. u. PROS. - 1 0.9000001*04 4.92570 l. 000000000 

DESIGN POINT COORD., SENSITIVITY P ACTORS, STD. D1V. PROM MEAN VALUES 
Z0, PROS, • 0.9000001*04 1.000000000 

WOD 0 . 5397(3 11*07 0.999997 -4.924923 

ALFA 0.99952731-05 -0.015057 0.102775 


CDF RESULTS 

z U PROBABILITY ITER. NO. 

0.900000001*04 0.692549541*01 0.100000001*01 0 


87 



gfun.mo-v 


Sun Jun 29 03*16*45 1997 


1 


CONDITION • NODS COMPONENT LAYER 


INCR 0 44 3 I 


I MEDIAN, KEAN, STD. DTV. 

0.6782718*04 0.6782711*04 0.2127548*03 

NUMBER OF PROBABILITY LEVELS 1 

LVL, Z0, u, PROB. ■ 1 0.800000E*04 5.72158 0.999559955 

DESIGN POINT COORD., SENSITIVITY FACTORS, STD. DEV. FROM KEAN VALUES 
Z0, PROB. • 0.0000001*04 0.999*55895 

■MOD 0.72744078*07 0.955886 -5.720537 

ALFA 0.96640208-05 -0.015088 0.086326 


CDF RESULTS 

I Q PROBABILITY ITER. NO. 

0.100000001*04 0.572151818*01 0.999999591*00 0 


gfun.mov Sun Jun 29 03*11*36 1997 1 


CONDITION I NODE COMPONENT LAYER 


INCR 0 44 3 1 


Z MEDIAN, MEAN, STD. DEV. 

0.7016048*04 0.7016048*04 0.2129948*03 

NUMBER OF PROBABILITY LEVELS 1 

LVL, Z0, u, PROB. ■ 1 0.8000008*04 4.61965 0.959998076 

DESIGN POINT COORD., SENSITIVITY FACTORS, STD. DEV. FROM MIAN VALUES 
Z0, PROB. • 0 . 800000E+04 0.99999B076 

EMOD 0 . 91474926*07 0.999B86 -4.619122 

ALFA 0.96322828-05 -0.015071 0.069622 


CDF RESULTS 

z U PROBABILITY ITER. NO. 

0.800000008*04 0.461964718*01 0.999998088*00 0 


gfun.mcrv Sun Jun 29 03*13*50 1997 1 


CONDITION I NODE COMPONENT LAYER 


INCR 0 44 3 1 


Z MEDIAN. MEAN. STD. DEV. 

0.7482711*04 0.7482718*04 0.2134541*03 

NUMBER OF PROBABILITY LEVELS 1 

LVL, Z0, u, PROB. • 1 0.8000008*04 2.42343 0.992312628 

DESIGN POINT COORD., SENSITIVITY FACTORS, STD. DEV. PROM MEAN VALUES 
Z0, PROB. - 0.8000008*04 0.992312(28 

EMOD 0.12880648*08 0.999887 -2.423155 

ALFA 0.95692448-05 -0.015038 0.036444 


CDF RESULTS 

l U PROBABILITY IT8R. NO. 

0 . 800000001*04 0.242342928*01 0.992312638*00 0 
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gfun.nov 


Tub Jul 1 15 i 39 s 56 1997 


1 


CONDITION f HODS COMPONENT LAYER 


I34CR 0 44 3 1 


2 MEDIAN. KEAN, STD. DEV. 

0 . 7714042*04 0 . 7716042+04 0.2136941*03 

NUKBW OF PROBABILITY LEVELS 1 

LVL, 20. U, PROS. • 1 0.9000002*04 1.32910 0.909043193 

D RSI 31 POINT COORD.. SENSITIVITY FACTORS, STD. DEV. FROM MEAN VALORS 
20, FROi. - 0.8000002*04 0.909043143 

WOO 0.14741292*01 0.999997 >1.329451 

ALFA 0.95379252-05 -0.015021 0.019941 


CDF RESULTS 
S 0 

0.900000001*04 0.132990122*01 


PROBABILITY 

0.909043142*00 


gfun.mov 


Tub Jul 1 15 : 36:47 1997 1 


CONDITION 


NODI COMPONENT LAYER 


/ 


incr 


o 


44 


1 


Z MEDIAN. KEAN, STD. DEV. 

0.7949382*04 0.7949382*04 0.2139242*03 

NUMBER OF PROBABILITY LEVELS 1 

LVL. 20, u. PROB. « 1 0.9000002*04 0.23444 0.593531990 

DESIGN POINT COORD., SENSITIVITY FACTORS, STD. DRV. FROM KKAN VALUES 
20, PROB. - 0.800000**04 0.593531990 

EMOD 0.1659776**08 0.999898 >0.234414 

ALFA 0.9504724E>05 >0.014959 0.003540 


CDF RESULTS 
2 U 

0.800000002*04 0.234440192*00 


PROBABILITY 

0.593531992*00 


gfun.BOT 


Tub Jul 1 16 : 17:30 1997 1 


CONDITION 


MODE COMPONENT LAYER 


INCR 0 44 


3 


1 


2 MEDIAN, KEAN, STD. DEV. 

0.7996041*04 0.7996042*04 0.2139642*03 

NUMBER OF PROBABILITY LEVELS 1 

LVL, 20, u, PROB. • 1 0.8000002*04 0.01849 0.507377547 

DE5IQ4 POINT COORD., SENSITIVITY FACTORS, STD. DEV. FROM KEAN VALUES 
20. PROB. * 0.1000002*04 0.507377567 

EMOD 0. 16968562+08 0.999887 -0.018492 

ALFA 0.95005272-05 -0.015003 0.000277 


2 

0. 800000001+04 


CDF RESULTS 

U PROBABILITY 

0.114937582-01 0.507377572*00 


ITER. 140. 


0 
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Appendix G 


DATA AND OUTPUT FILES FOR CALCULATION OF UNION 


PROBABILITY 


(Chapter IV) 


temp.d&t Tue Jul 1 16:23:17 1997 

7, 2, 1 

X, 17000000.0,1700000.0 
2, 0.0000095,0.0000019 

1, 6.825696 

2, 5.721588 

3, 4.619647 

4, 2.423429 

5, 1.328801 

6, 0.236641 

7, 0.018497 

1, 6.825696 0.999887 -0.015057 

2, 5.721588 0.999886 -0.015088 

3, 4.619647 0.999886 -0.015071 

4, 2.423429 0.999887 -0.015038 

5, 1.328801 0.999887 -0.015021 

6, 0.236641 0.999888 -0.014959 

7, 0.018497 0.999887 -0.015003 


***** BOUNDS ***** 

Lower bound * 0.4926211 
Upper bound = 0.4926211 


90 



BIBLIOGRAPHY 


[1] Mahadevan, S., Mehta, S., Tryon, R. G. and Cruse, T. A., “System Reliability 
Design Analysis of Propulsion Structures”, Presented at the International Gas 
Turbine and Aeroengine Congress and Exposition, Cincinnati, Ohio, May, 1993. 

[2] Ang, A. H-S. and Tang, W. H., “Probability Concepts in Engineering Planning 
and Design, Vol.I - Basic Principles”, John Wiley & Sons, Inc., 1975. 

[3] Ang, A. H-S. and Tang, W. H., “Probability Concepts in Engineering Planning 
and Design, Vol.II - Decision, Risk, and Reliability”, John Wiley & Sons, Inc., 
1984. 

[4] Pages, A. and Gondran, M., “System Reliability: Evaluation and Prediction in 
Engineering”, Springer - Verlag, 1986. 

[5] Tung, A. T. Y. and Kiremidjian, A. S., “Application of System Reliability The- 
ory in the Seismic Analysis of Structures”, Earthquake Spectra, Vol.8, No.3, 
1992. 

[6] Incropera, F. P. and DeWitt, D. P., “Introduction to Heat Transfer, 2nd ed.”, 
John Wiley Sons, Inc., 1990. 

[7] Thomas, L. C., “Heat Transfer”, Prentice-Hall, Inc., 1992. 

[8] Logan, D. L., “A First Course in the Finite Element Method”, PWS Publishers, 
1986. 

[9] Kraus, A. D. and Bar-Cohen, A., “Thermal Analysis and Control of Electronic 
Equipment”, Hemisphere Publishing Corporation, 1983. 

[10] Janna, W. S., “Introduction to Fluid Mechanics”, PWS Publishers, 1983. 

[11] Southwest Research Institute, “NESSUS Reference Manual, Version 1.0”, July, 
1991. 

[12, Southwest Research Institute, “A Short Course on Probabilistic Structural Anal- 
ysis Methods & NESSUS Workshop, Vol. II”, 1990. 


91 



[13] Dias, J. B., “The NESSUS Finite Element Code - User’s Manual”, Southwest 
Research Institute, 1989. 


92 



